AD-A119  9?9  NAVAL  CIVIL  ENGINEERING  LAB  PORT  HUENEME  CA  F/G  13/10 

STATNOOR  -  A  SINGLE-POINT  MOORING  STATIC  ANALYSIS  PROGRAM. (U) 

,  JUN  82  J  V  COX 

UNCLASSIFIED  NCEL-TN-1634  NL 


DTK  FILE  COPY  ADAH"79 


TN  NO:  N-1634 


jlj.  p.  STATMOOR  -  A  SINGLE-POINT  MOORING 
I  I  I  UC.  STATIC  ANALYSIS  PROGRAM 


AUTHOR:  j.  v.  cox 


DATE:  June  1982 


SPONSOR:  Naval  Facilities  Engineering  Command 


PROGRAM  NO:  yfs?. 556 .091.01.403 


NAVAL  CIVIL  ENGINEERING  LABORATORY 
PORT  HUENEME,  CALIFORNIA  93043 

Approved  for  public  release;  distribution  unlimited. 


82  10  0 7 


w 

12 


METRIC  CONVERSION  FACTORS 


’tin  -  2-54  Fof  othar  exact  conwrtionc  and  mor*  detailed  tablet,  mm  NBS 

Mite.  Put*.  286.  Unit*  of  Wtight*  and  Maaturac,  Prica  $2.25,  SD  Catalog  No  Cl  3.10:286. 


I 

Unclassified 


SECURITY  CLASSIFICATION  OF  TwiS  PAGE  r*hen  Dote  Fniw.ed) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  REPORT  NUMBER  2 

TN-1634  DN78701 1 

J  RECIPIENT'S  CATALOG  NUMBER 

4  title  (mod  Subtitle) 

STATMOOR  -  A  SINGLE-POINT  MOORING 

STATIC  ANALYSIS  PROGRAM 

S  TYPE  OF  REPORT  A  PERIOD  COVERED 

Not  final;  Aug  1981  —  Apr  1982 

6  PERFORMING  ORG  REPORT  NUMBER 

7  AuTHORfa, 

J.  V.  Cox 

6  CONTRACT  OR  GRANT  NuMBCRfai 

9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

NAVAL  CIVIL  ENGINEERING  LABORATORY 

Port  Hueneme,  California  93043 

10  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  6  WORK  UNIT  NUMBERS 

62759N; 

YF59. 5  56.091. 01. 403 

J>  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Facilities  Engineering  Command 

Alexandria,  Virginia  22332 

»2  REPORT  DATE 

June  1982 

1)  NUMBER  OF  PAGES 

39 

14  MONITORING  AGENCY  NAME  *  ADDRESSfrf  different  front  Controlling  Office) 

is  SECURITY  CLASS  (ol  thie  report) 

Unclassified 

»S  DISTRIBUTION  STATEMENT  (ol  thin  Report) 

Approved  for  public  release;  distribution  unlimited. 

»7  DISTRIBUTION  STATEMENT  (of  the  obetrect  entered  in  Block  70.  it  diUetont  tram  Rep ort) 


If.  SUPPLEMENTARY  NOTES 


If.  KEY  WORDS  (Contlnuo  on  rererte  aide  it  neceeamry  ond  Identity  by  block  number) 

Moorings,  single-point  moorings,  wind  forces,  current  forces,  mooring  analysis  computer 
program,  hawser  loads. 


20  ABSTRACT  (Continue  on  rovotmo  aide  It  nacaaaoty  ond  identity  by  block  number) 

/STATMOOR  is  a  sutic  mooring  analysis  program  written  in  BASIC  language  and  is 
one  program  in  a  hierarchy  of  programs  developed  at  the  Naval  Civil  Engineering  Laboratory 
for  mooring  analysis.  STATMOOR  analyzes  the  static  response  of  a  single-point  moored 
vessel  and  hawser.  The  MENU  arrangement  of  the  program  lends  itself  to  a  user-oriented 
conversational  mode.  The  user  has  the  option  to  enter,  review,  edit  input,  and  obtain 

(continued) 

DO  ,'23",  1473  eo.T.oNOF  .Nov.j.so.soLETt  Unclassified 

SECURITY  CLASSIFICATION  OP  TM1*  PACE  (•»!•#»  DmttTEntOfOd) 


_ Unclassified _ 

ItemiTy  CLASSIFICATION  OF  THIS  PAGtmtn  Dm  Bn ItnJI 

20.  Continued 

calculated  results  in  printed  tabular,  video  tabular,  or  video  graphics  form.  Steady  current, 
wind,  and  wave  loads  are  considered.  Wind  load  estimates  are  considered  to  be  as  accurate 
as  the  user’s  knowledge  of  the  wind  environment;  current  and  wave  loads  are  in  a  pre¬ 
liminary  form  and  merit  further  refinement.  STATMOOR  was  written  to  demonstrate 
the  utility  and  ease  of  use  of  conversational  mode  programs  and  the  potential  for  computer 
programs  to  replace  bound  design  manuals.  It  is  unique  in  that  it  incorporates  recently 
developed  wind  load  information,  can  confidently  be  used  with  little  or  no  training,  and 
is  easily  adaptable  to  most  desk  top  computers. 


Library  Card 

Naval  Civil  Engineering  Laboratory 

STATMOOR  -  A  SINGLE-POINT  MOORING  STATIC 

ANALYSIS  PROGRAM,  by  J.  V.  Cox 

TN-1634  39  pp  illus  June  1982  Unclassified 

1.  Current  forces  2.  Conversational  mode  program  I.  YF59.556.091.01.403 

STATMOOR  is  a  static  mooring  analysis  program  written  in  BASIC  language  and  is  one 
propam  in  a  hierarchy  of  programs  developed  at  the  Naval  Civil  Engineering  I  aboratory  for 
mooring  analysis.  STATMOOR  analyzes  the  static  response  of  a  single-point  moored  vessel  and 
hawser.  The  MENU  arrangement  of  the  program  lends  itself  to  a  user-oriented  conversational 
mode.  The  user  has  the  option  to  enter,  review,  edit  input,  and  obtain  calculated  results  in 
printed  tabular,  video  tabular,  or  video  graphics  form.  Steady  current,  wind,  and  wave  loads 
are  considered.  Wind  load  estimates  are  considered  to  be  as  accurate  as  the  user's  knowledge  of 
the  wind  environment;  current  and  wave  loads  are  in  a  preliminary  form  and  merit  further 
refinement.  STATMOOR  was  written  to  demonstrate  the  utility  and  ease  of  use  of  conversa¬ 
tional  mode  programs  and  the  potential  for  computer  programs  to  replace  bound  design 
manuals.  It  is  unique  in  that  it  incorporates  recently  developed  wind  load  information,  can 
confidently  be  used  with  little  or  no  training,  and  is  easily  adaptable  to  most  desk  top  computers. 


SCCUKirr  CLASSIFICATION  or  THIS  FAS erNAm  OM«  FMI.HI 


u*4i 


CONTENTS  i 

| 

! 

Page 


INTRODUCTION  .  1 

i 

OPERATION  .  1 

Master  Menu  Options  .  1 

Input  . 

Angle  Conventions  and  Coordinate  System 
Computations  . 

BENEFITS  AND  EXTENSIONS  .  3 

REFERENCES  .  4 

APPENDIXES 

A  -  Equations  .  6 

B  -  Block  Flow  Chart,  Sample  Output,  and  Listings  .  16 


CN  CM  CO 


INTRODUCTION 


The  Naval  Facilities  Engineering  Command’s  (NAVFAC)  mission  includes 
the  design  and  construction  of  fixed  ocean  facilities.  To  support  this 
aspect  of  the  mission,  the  Naval  Civil  Engineering  Laboratory  (NCEL)  is 
developing  a  hierarchy  of  mooring  analysis  programs,  which  include  both 
static  and  dynamic  analysis.  This  report  describes  the  use  and  contents 
of  STATMOOR,  which  analyzes  the  static  response  of  a  single  point  moored 
vessel  and  hawser.  It  calculates  the  hawser  load,  hawser  angle,  and 
ship  angle  for  steady  current,  wave,  and  wind  loads.  Wind  load  estimates 
are  considered  to  be  accurate;  current  and  wave  loads  are  in  a  preliminary 
form  and  merit  further  refinement. 

The  program  is  written  in  a  conversational  mode  which  is  described 
in  the  operating  section.  The  benefits  and  extensions  of  this  work 
regarding  Navy  mooring  applications  are  discussed  later.  The  equations 
used  in  STATMOOR  for  the  current,  wave,  and  wind  loads  on  the  ship  are 
in  Appendix  A.  Appendix  B  contains  a  program  listing,  operation  flow 
chart,  and  example  printout. 


OPERATION 

The  information  in  this  section  is  overviewed  in  the  explanation 
page  of  the  program  (Appendix  B) ,  but  is  described  in  more  detail  here. 

Master  Menu  Options 

The  program  is  menu  driven  to  lend  to  its  user  oriented  conver¬ 
sational  mode.  The  master  menu  allows  the  user  to  access  the  various 
program  options.  When  the  following  master  menu  is  presented  the  user 
is  prompted  to  enter  the  number  of  the  option  choice.  Options  1  through 
7  return  control  to  the  master  menu  after  execution  (as  diagrammed  in 
Appendix  B) . 

1.  Input  vessel  characteristics. 

2.  Input  environmental  conditions. 

3.  Review/edit  vessel  characteristics. 

A.  Review/edit  environmental  conditions. 

5.  Calculate. 

6.  Tabular  output. 

7.  Graphic  display. 

8.  Exit  program. 

Enter  Number  of  above  choice? 


1.  Input  Vessel  Characteristics  -  Allows  the  user  to  enter  information 
about  the  vessel  by  presenting  the  user  vith  a  series  of  questions/choices. 
The  answers  provide  the  necessary  input. 


2.  Input  Environmental  Conditions  (Similar  to  Option  1)  -  Allows  the 
user  to  enter  information  about  the  environmental  conditions  (current, 
wave,  and  wind). 

3.  Review/Edit  Vessel  Characteristics  -  Allows  the  user  to  review  and 
change  any  of  the  vessel  characteristics. 

4.  Review/Edit  Environmental  Conditions  -  Allows  the  user  to  review  and 
change  any  number  of  the  environmental  conditions.  [NOTE:  Options  3 
and  4  allow  parametric  studies  without  redundant  input.] 

5.  Calculate  -  Performs  the  calculations  of  the  forces  on  the  ship  and 
the  equilibrium  conditions.  Some  loading  conditions  yield  two  stable 
solutions,  so  the  user  has  the  option  to  search  for  multiple  stable 
solutions.  Another  user  option  is  the  ability  of  the  program  to  determine 
the  direction  of  current  loading  (within  10  degrees)  which  causes  the 
maximum  hawser  load.  A  brief  summary  of  the  output  can  either  be  printed 
or  displayed  on  the  video. 

6.  Tabular  Output  -  Allows  the  user  to  obtain  either  printed  or  video 
output.  Video  output  consists  of  a  table  of  force  coefficients,  a  taole 
of  the  individual  environmental  forces  on  the  vessel,  and  equilibrium 
conditions  (hawser  load,  hawser  angle,  and  ship  angle).  Printed  output 
consists  of  all  the  video  output  plus  all  the  vessel  and  environmental 
input  (see  Apnendix  B) . 

7.  Graphic  Display  -  Gives  a  graphic  display  of  the  incident  current, 
wave,  and  wind  directions  as  well  as  equilibrium  angle  of  the  ship  and 
the  hawser. 

8.  Exit  Program  -  Exits  the  program  and  returns  control  to  the  operating 
system. 

Input 


The  input  consists  of  numerical  values,  letters,  and  descriptions, 
pi  'mpted  by  the  program: 

•  Numerical  values  must  be  typed  with  no  alpha  characters  (e.g. , 
no  letters) . 

•  Letters  are  used  as  input  to  subjective  questions  and  to 
control  the  flow  of  the  program.  At  several  points  the  program 
asks  a  question  followed  by  the  characters  (1/N),  which  means 
the  user  should  enter  Y  for  yes  or  N  for  no. 

•  Descriptions  are  entered  in  the  input  sections  of  the  program 
to  document  the  problem  being  solved.  Alpha  and  numeric 
characters  can  be  used  in  the  description. 

All  input  must  be  followed  by  the  return  key. 

Angle  Conventions  and  Coordinate  System 

There  are  two  angle  conventions  in  the  program,  absolute  and  local. 
The  fixed,  absolute  angle  is  zero  at  twelve  o'clock  and  increases  counter¬ 
clockwise.  It  is  used  to  enter  the  incident  angle  of  the  environmental 
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loads  and  output  the  equilibrium  angles  of  the  ship  and  hawser.  The 
local  angle  rotates  relative  to  the  vessel  direction.  Head  loading  is 
zero  degrees  and  the  angle  increases  counterclockwise.  The  local  angle 
is  ured  for  calculation  of  loads  and  output  in  the  force  table.  The 
local  coordinate  system  is  also  relative  to  the  vessel.  The  positive 
X-axis  is  at  the  local  angle  of  180  degrees,  toward  the  stern;  positive 
Y-axis  is  at  the  local  angle  of  270  degrees,  toward  the  starboard.  The 
local  angle  convention  and  coordinate  system  are  shown  in  Figure  1. 


local  angle  for  current, 
wave,  or  wind  attack 


90 


y  270 


180 


Figure  1.  Local  angle  and  coordinate  system. 


Computations 

STATMOOR  considers  steady  current,  wind,  and  wave  loads.  Each  of 
these  environmental  loads  is  broken  into  a  longitudinal  force,  lateral 
force,  and  yaw  moment.  Wind  load  computations  consider  the  wind  gradient 
as  well  as  the  vessel  geometry  (vessel  group).  The  wind  loads  are  based 
on  recommendations  from  Reference  1.  The  current  loads  include  form, 
friction,  and  propeller  drag,  and  are  based  on  References  2  and  3.  The 
wave  drift  loads  consider  the  wave  spectrum  and  drift  force  operators 
and  are  based  on  References  2  and  4  through  8.  The  equations  are  presented 
in  Appendix  A. 


BENEFITS  AND  EXTENSIONS 

STATMOOR  was  written  to  demonstrate  the  feasibility  of  conver¬ 
sational  mode  programs  and  is  one  of  a  hierachy  of  programs  developed 
for  mooring  analysis  at  NCEL.  Its  simplicity  and  ease  of  input  allow  it 
to  be  used  with  little  or  no  instructions. 

STATMOOR  is  unique  for  the  following  reasons: 

(1)  First  program  to  incorporate  recently  developed  wind  load 
information  (Ref  1). 
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(2)  Conversational  to  the  extent  it  can  be  confidently  used  with 
little  or  no  training. 

(3)  Easily  adaptable  to  most  desk  top  computer. 

The  use  of  this  computer  program  has  advantages  over  both  hand 
calculations  and  tabulated  results,  as  outlined  below: 

(1)  Reduced  calculation  time  (especially  to  determine  worst  case 
loadings) . 

(2)  Less  chance  of  numerical  error. 

(3)  No  chance  of  procedural  error. 

(4)  Ease  of  parametric  studies. 

(5)  Ability  to  quickly  check  the  effect  of  input  uncertainties  on 
the  results. 

(6)  More  variables  can  be  considered  than  are  practical  for  tabulated 
results . 

(7)  Advances  in  state-of-the-art  can  be  incorporated  with  ease. 

In  its  present  form  the  output  equilibrium  condition  can  be  used  to 
enter  more  sophisticated  computer  programs  which  do  not  easily  converge 
on  the  static  solution  or  to  enter  programs  which  determine  the  response 
of  the  complete  mooring. 

A  conversational  program  could  be  expanded  to  include  design  aids, 
such  as  ultimate  strength  as  a  function  of  cable  diameter  relationships 
(Ref  9).  The  result  would  be  a  very  cost-effective  tool  for  preliminary 
design  and  cost  estimating.  With  further  extensions  this  type  of  program 
could  become  a  computer-based  mooring  design  manual,  as  an  alternative  to 
the  present  DM  26  (Ref  10) . 


REFERENCES 

1.  Naval  Civil  Engineering  Laboratory.  Technical  Memorandum  M-44-81-7: 
Wind  induced  steady  loads  on  moored  ships,  by  R.  Owens  and  P.  A.  Palo. 
Port  Hueneme,  Calif.,  1981. 

2.  G.  F.  M.  Remery  and  G.  Von  Oortmerssen.  "The  mean  wave,  wind,  and 
current  forces  on  offshore  structures  and  their  role  in  the  design  of 
mooring  systems,"  paper  presented  at  Offshore  Technology  Conference, 
Houston,  Tex.,  1973.  (Paper  no.  OTC  1741) 

3.  Hydronautics ,  Inc.  Technical  Report  7096-1:  Forces  on  ships  moored 
in  protected  waters,  by  R.  Altman.  Laurel,  Md.,  1971. 

4.  0.  M.  Faltinsen  and  A.  E.  Loken.  "Drift  forces  and  slowly  varying 
forces  on  ships  and  offshore  structures  in  waves,"  Norwegian  Maritime 
Research,  vol  6,  no.  1,  1978. 


4 


5.  Netherlands  Ship  Model  Basin.  Wave  drifting  forces,  by  J.  A.  Pinkster. 
1979. 


6.  S.  K.  Chakrabarti.  "Steady  and  oscillating  drift  forces  on  floating 
objects,"  ASCE  Journal  of  the  Waterway,  Port,  Coastal  and  Ocean  Division, 
vol.  106,  no.  WW2,  May  1980. 

7.  0.  M.  Faltinsen  and  F.  C.  Michelsen.  "Motions  of  large  structures 

in  waves  at  zero  Froude  number,"  paper  presented  at  International  Symposium 
on  the  Dynamics  of  Marine  Vehicles  and  Structures  in  Waves,  University 
College,  London,  Apr  1974. 

8.  Stevens  Institute  of  Technology,  Davidson  Laboratory.  Feasibility 
of  conducting  current  load  tests  with  ship  models.  Hoboken,  N.  J., 

9  Jan  1981.  (Informal  Document) 


9.  Hydrospace  Research  Corporation.  Technical  Report  270:  Designer's 
guide  for  deep-ocean  ship  moorings.  Rockville,  Md.,  1970. 

10.  Naval  Facilities  Engineering  Command.  Design  Manual  DM  26.6: 
Mooring  design  physical  and  experimental  data,  Preliminary  draft. 
Alexandria,  Va.,  1981. 


11.  Naval  Civil  Engineering  Laboratory.  Technical  Note:  An  assessment 
of  state-of-the-art  methods  for  calculating  current  loads  on  moored 
ships,  by  Paul  Palo  and  Robert  Owens.  Port  Hueneme,  Calif.,  1982. 


12.  B.  Molin.  "Computations  of  wave  drift  forces,"  paper  presented  at 
Offshore  Technology  Conference,  Houston,  Tex.,  1979.  (Paper  no.  0TC  3627) 


13.  J.  C.  Hungerford.  "Graphic  manipulations  using  matrices,"  Byte, 
Sep  1978. 


5 


stem . * 

— * _  _ 


Appendix  A 
EQUATIONS 


This  section  of  the  report  reviews  the  equations  used  in  STATMOOR 
for  the  current,  wave,  and  wind  loads  on  the  ship.  The  various  experi¬ 
ments  and  analytical  methods  examined  for  this  report  produced  a  large 
scatter  in  the  forces  induced  on  a  moored  ship  by  currents  and  waves. 
The  user  must  keep  this  uncertainty  in  current  and  wave  forces  in  mind. 
The  wind  forces  are  based  on  the  recommendations  in  Reference  1  and  are 
as  accurate  as  the  user's  knowledge  of  the  actual  wind  environment. 


CURRENT  LOADS  ON  VESSELS 


Currents  produce  longitudinal  (F  )  and  lateral  forces  (F  )  as 
well  as  a  yaw  moment  (M  )  about  the  center  of  gravity  of  the  vessel. 
Figure  A-l  shows  the  scatter  obtained  for  lateral  current  loads  by 
various  methods  in  an  NCEL  study  (Ref  11). 

Longitudinal  Current  Loads 

The  longitudinal  force  consists  of  forces  due  to  form,  friction, 
and  propeller  drag  as  follows: 

F  =  F  ,  +  F  ,  .  +  F  (A-l) 

xc  x  form  x  friction  x  prop 

Pressure  drag,  commonly  called  form  drag,  is  characterized  by 
separation  and  the  formation  of  a  downstream  wake.  Form  drag  is  given 
by  the  following  equation: 


F  ,  =  1/2  p  V2  C  .  A,  cos30  (A 

x  form  w  c  xcb  d  c 

where  Pw  =  mass  density  of  water  (1.9876  slug/ft3) 

=  average  current  speed  (ft/ sec) 

C  .  =  longitudinal  current  form  drag  coefficient 
xc  =0.1  (NCEL  estimate) 

A.  =  end-projected  area  below  the  waterline  of  the  vessel 
=  Beam  x  Draft  (ft2) 

6c  =  local  angle  of  the  current 


Skin  friction  drag  is  the  resultant  of  tangential  shear  stresses 
due  to  viscosity  and  velocity  gradients  at  the  boundary  surface.  Friction 
drag  is  given  by  the  following  equation  (Ref  2): 


F 

x 


friction 


1/2  p  V2  C  S  cos30 
w  c  xca  c 


(A-3) 


where 


C  =  longitudinal  skin  friction  coefficient 
XCa  =  0.075/(Log  Rn  -  2)2 

Rn  =  Reynolds  Number  =  L  cos8/\> 

L  =  waterline  length  of  the  vessel  (ft) 

v  =  kinematic  viscosity  of  water  (1.4  x  10  3  ft2/sec) 

S  =  wetted  surface  area  (ft2) 

=  (1.7  TL)  +  (35  D/T) 

T  =  draft  of  vessel  (ft) 

D  =  displacement  of  ship  (long  tons) 


Propeller 
locked  shaft 


drag  is  the  form  drag  of  the  vessel's  propeller  with  a 
Propeller  drag  is  given  by  the  following  equation  (Ref  4): 


F 

x  prop 


1/2  p  V2 
rw  c 


A 

P 


C 

prop 


cos 


(A-4) 


where  =  propeller  expanded  (or  developed)  blade  area  (ft.2) 

Cprop  =  ProPe^^er  dtag  coefficient  (assumed  to  be  1) 

Also  from  Reference  3: 


A  total  projected  propeller  area  , A  c  A 

p  "  1.067  -  0.229  p/d  1  ^ 

where  p/d  -  propeller  pitch  to  diameter  ratio  (assumed  to  be  1) 

The  angle  function  for  the  form,  friction,  and  propeller  drag  were 
taken  from  apparent  trends  in  experimental  data  from  experiments  conducted 
at  the  Stevens  Institute  of  Technology,  Hoboken,  N.J.  (Ref  8). 

To  limit  the  amount  of  detailed  input  to  the  program,  an  attempt 
was  made  to  express  the  total  projected  propeller  area  for  various 
vessels  in  terms  of  known  quantities.  Figure  A-2  shows  the  area  ratio 
(length  x  beam/total  projected  propeller  area)  versus  length  x  beam  for 
the  program's  five  major  vessel  groups.  Each  data  point  represents  a 
class  of  ship.  The  area  ratio  value  given  to  represent  each  group  of 
data  points  considers  the  number  of  ships  in  each  class.  Cruisers, 
cargo,  and  tankers  are  represented  by  a  single  area  ratio.  Destroyers 
and  carriers  each  are  represented  by  two  distinct  groups  as  a  function 
of  L  x  B.  The  total  projected  propeller  area  can  be  given  in  terms  of 
the  area  ratio  as  follows: 


Total  projected  propeller  area 


L  x  B 

area  ratio 


(A-6) 
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The  simplicity  in  this  method  of  determining  propeller  drag  is 
justified  at  this  time  because  of  the  overall  uncertainty  in  the  current 
load  equations.  Its  arrangement  in  the  program  lends  itself  to  later 
refinement . 

Lateral  Current  Loads 


The  lateral  current  load  on  the  vessel  is  due  primarily  to  form 
drag.  The  lateral  current  load  is  calculated  by  the  following  equation 
(Ref  2): 

F  =  1/2  p  V2  C  A  K  (A-7) 

yc  rc  c  yc  c 

where  A  =  side-projected  area  below  the  waterline  of  the  vessel 
C  =  length  x  draft  (ft2) 

C  =  lateral  current  load  coefficient 
yc 

K  =  depth  to  draft  factor 

The  lateral  current  load  coefficient  was  obtained  from  one  source, 
based  on  tanker  model  tests.  It  is  presented  as  a  Fourier  series  (Ref  2) 

5 

C  =  1  b  sin(n0  )  (A-8) 

yc  ,  n  c 

n=l 

in  which  b,  =  0.908  b.  =  0 

1  4 

b  =  0  b  =  -0.033 

b3  =  -0.116 

The  depth  to  draft  factor  in  Equation  A-7  compensates  for  shallow 
water  effects.  The  following  equation  was  fitted  to  the  curve  in 
Reference  2: 

K  =  1  +  - *-= —  (A-9) 

(WD/T)  -1 

where  WD  =  water  depth  (ft) 

T  =  vessel  draft  (ft) 

Yaw  Moment  Due  to  Currents 


The  yaw  moment  on  the  vessel  due  to  currents  is  calculated  by  the 
following  equation  (Ref  2): 


M 

c 


1/2  p  V2 
rc  c 


C 

xyc 


A 


c 


L  K 


where  C  =  current  yaw  moment  coefficient 
xyc  3 

K  =  depth  to  draft  factor  (Equation  A-9) 


(A- 10) 
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The  current  yaw  moment  coefficient  was  obtained  from  the  same 
source  as  the  lateral  load  coefficient.  Expanded  in  Fourier  series  the 
coefficient  is  given  as  follows: 


5 

C  =  S  bn  sin  n8  (A— 11) 

xyc  n=l 

in  which  b^  =  -0.0252 
b2  =  -0.0904 
b3  =  0.0032 

WAVE  DRIFT  LOADS  ON  VESSELS 

The  wave  drift  force  is  a  result  of  the  vessel  acting  as  an  obstruc¬ 
tion  to  the  transport  of  momentum  in  the  wave  field.  The  method  used  to 
determine  the  steady  drift  forces  (longitudinal,  lateral,  and  yaw  moment) 
assumes  the  drift  force  operators  on  the  vessel  as  a  function  of  the 
wave  frequency  is  known.  These  operators  (transfer  functions)  were 
based  primarily  on  model  tests.  The  transfer  functions  are  plotted 
versus  a  nondimens ional  wavelength.  The  program  converts  each  function 
to  the  angular  frequency  domain  by  the  following  deep  water  relationship: 


b^  =  0.0109 

b5  =  o.oon 


u> 


(A-12) 


There  is  a  larger  scatter  in  the  available  data  and  the  effect  of 
water  depth  to  draft  is  not  considered.  Figure  A-3  shows  lateral  drift 
forces  on  a  tanker  for  headings  of  90  and  45  degrees  (Ref  12).  From  this 
one  model  test  it  appears  the  lateral  drift  force  varies  with  the  sine 
squared  of  the  incident  angle.  In  a  similar  manner  the  longitudinal 
drift  force  is  assumed  to  vary  with  cosine  squared  of  the  incident 
angle. 

Wave  Spectra 

Two  wave  spectra  are  available  in  the  program.  Both  are  of  the 
following  form: 


S(u>)  =  A/u)5  e_(B/u,4)  (A-13) 

where  u>  is  the  wave  freqency. 

If  the  user  indicates  the  wave  period  is  known  (as  in  fetch  limited 
anchorages)  then  the  Bretschneider  wind  wave  model  is  used,  where: 

A  =  262.5  Hz/T4 
B  =  1050/T4 

H  =  significant  wave  height  (ft) 

T  =  significant  wave  period  (sec) 
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If  the  user  indicates  the  wave  period  is  not  known  then  the  Pierson- 
Moskowitz  wind  wave  model  is  used,  where: 


B  =  33.56/H2 

Longitudinal  Wave  Drift  Loads 


Figure  A-4  shows  representative  nondimensional  drift  force  coeffi¬ 
cients  versus  the  nondimensional  wavelength  for  longitudinal  loads. 

There  is  a  large  scatter  in  the  drift  force  values.  Additional  parameters 
such  as  the  bow  angle  should  possibly  be  included  in  the  nondimensional 
plot.  The  solid  curve  represents  the  transfer  function  used  in  the 
program,  where  as  the  dotted  and  dashed  curves  represent  experimental 
and  theoretical  values  respectively.  The  longitudinal  drift  force  is 
obtained  with  the  following  equation  (Ref  2): 


(  /  Fxd(u°  \ 

Fxd  =  Pw«  L  J  S(w>  — - - 2 

o  \1/2  pu>8  •  / 


■jldu)  cos  0d 


(A-14) 


where  g  =  acceleration  of  gravity  (32.2  ft/sec2) 

Fxd(w)/(l/2  g  a2)  =  longitudinal  drift  force  operator 
8d  =  local  angle  of  the  waves 
a  =  wave  amplitude 

Lateral  Wave  Drift  Loads 

Figure  A-5  shows  representative  nondimensional  drift  force  coeffi¬ 
cients  for  lateral  loads.  The  scatter  is  not  as  pronounced  as  it  was 
for  longitudinal  coefficients.  The  transfer  function  in  STATM0OR 
(represented  by  a  solid  line)  is  a  five-line  sequence  representation  of 
the  drift  forces  on  a  flat  plate.  The  lateral  drift  force  is  obtained 
by  the  following  equation  (Ref  2): 


i  v/2  V  ■  / 


Jdu)  sin3  ed 


(A-15) 


where  F  .(w)/(l/2  p  g  a2)  is  the  transfer  function  for  lateral  drift 
forces . 

NOTE:  Figures  A-4  and  A-5  have  different  scales  than  the  equations. 
Yaw  Moment  Due  to  Waves 

Figure  A-6  shows  the  nondimensional  yaw  drift  moment  for  a  heading 
of  45  degrees.  Only  two  theoretical  curves  (dashed)  are  shown  with  the 
program's  (solid)  transfer  function.  The  maximum  yaw  moment  due  to  wave 
drift  forces  was  assumed  to  be  at  45  degrees  and  periodic  over  180  degrees. 
The  resulting  yaw  moment  equation  is  given  by: 


i 


Md  =  Pu»8  1,2  /  S(w)| 


Md(ui) 


,1/2  Pu,8  L  a 


2  idu>  sin(20d) 


(A-16) 


where  Md(u»)/(l/2  p^g  L  a2)  is  the  transfer  function  for  yaw  moment  due 
to  waves . 


WIND  LOADS  ON  VESSELS 


All  methods  for  determining  vessel  wind  loads  were  taken  from 
Reference  1,  which  summarized  results  based  on  experimental  wind  load 
measurements  for  over  40  ships.  The  methods  in  the  report  do  not  take 
the  simplistic  approach  for  representing  the  forces  on  all  vessels  by 
three  force  versus  angle  curves.  The  force  curves  are  instead  a  function 
of  the  vessel  characteristics.  Among  some  of  the  items  considered  are: 

•  The  relative  areas  and  positions  of  the  hull  and  superstructure 
in  the  wind  gradient. 

•  The  nature  of  the  superstructure:  distinct  or  distributed. 

•  The  relative  longitudinal  position  of  the  superstructure. 

This  last  item  is  illustrated  in  Figure  A-7,  where  a  trend  in  the 
wind  yaw  moment  can  be  observed  as  the  location  of  the  superstructure 
changes . 

Longitudinal  Wind  Loads 

The  longitudinal  wind  load  is  a  strong  function  of  the  ship  type. 
Both  the  longitudinal  wind  drag  coefficient  and  the  normalized  shape 
function  (the  manner  in  which  the  load  changes  as  a  function  of  the 
angle)  are  based  upon  subjective  information  about  the  vessel.  The 
longitudinal  wind  load  is  given  by  the  equation: 


F  =  1/2  p.  V2  A  C  f  (6  )  (A-17 

xw  A  w  x  xw  x  w 

_  2 

where  p^  =  the  density  of  air  (2.4  x  10  slug/ft3) 

=  wind  velocity  at  a  reference  height  of  33  ft  (ft/sec) 

A  =  longitudinal  projected  vessel  area  above  the 
waterline  (ft2) 

=  longitudinal  wind  drag  coefficient 

f  (0  )  =  normalized  longitudinal  shape  function  dependent  on 
incident  wind  angle  (0^) 


For  longitudinal  loads  the  shape  function  can  be  one  of  two  forms. 
For  vessels  with  single  distinct  superstructures  and  for  hull  dominated 
ships  the  shape  function  takes  the  following  form: 


f  (8  )  =  cos  ♦ 

x  w 


(A-18) 
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Vessels  with  distributed  superstructures  have  a  shape  function  of 
the  following  form: 


f  (6  ) 

x  w 


sin  y  ”  (sin5y)/10 
1  -  1/10 


(A-19) 


In  both  shape  functions  the  angles  4>  and  y  are  a  function  of  the 
incident  wind  angle  and  the  longitudinal  location  of  the  superstructure. 
Equation  A-19  gives  a  flattened  out  shape  function  around  the  head  wind 
conditions . 


Lateral  Wind  Loads 


Calculation  of  the  lateral  wind  load  is  more  straightforward  than 
the  longitudinal  wind  load.  As  mentioned  in  the  previous  section,  the 
wind  velocity  used  in  the  wind  load  calculations  is  for  a  height  of 
33  feet  (10  meters).  The  user  enters  the  known  wind  velocity  and  cor¬ 
responding  height  at  which  it  is  measured,  and  the  program  determines 
the  calculation  wind  speed  from  the  wind  gradient  relationship: 


(A-20) 


where  =  reference  height  at  which  the  wind  velocity  was 

measured  (ft) 

VR  =  user-supplied  wind  velocity  at  the  reference  height  (knots) 

=  wind  velocity  at  33  feet  (knots,  program  converts  it  to 
ft/sec  for  calculations) 


The  lateral  wind  drag  coefficient  is  obtained  by  considering  the 
position  of  both  the  hull  and  the  superstructure  within  the  wind  gradient. 
The  lateral  wind  load  is  given  by  the  equation: 


F 

yw 


1/2  pA  V*  A  C  f  (6  ) 
A  w  y  y  y  w 


(A-21) 


where 


A 


y 

C 


y 


=  lateral  project  vessel  area  above  the  waterline  (ft2) 
=  lateral  wind  drag  coefficient,  where 


(A-22) 


where 


f 


y 


Ag  =  lateral  projected  superstructure  area  (ft2) 

A^  =  lateral  projected  hull  area  above  the  waterline  (ft2) 
Vg  =  average  wind  velocity  over  the  superstructure  (ft/sec) 
Vy  =  average  wind  velocity  over  the  hull  (ft/sec) 

(8^)  =  normalized  lateral  shape  function 
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Both  (V  /V^)  and  (V  /V^)  are  obtained  by  integrating  the  normalized 
wind  gradient  curve  (similar  to  Equation  A-20).  The  curve  is  integrated 
between  the  top  and  bottom  height  of  the  assumed  rectangular  superstructure 
or  hull  area. 

The  following  normalized  shape  function  was  used  for  all  vessels: 


f  (6  ) 

y  w 


sin  6  -  sin(56  )/20 

_ w _  w 

1  -  1/20 


(A-23) 


This  is  similar  to  the  longitudinal  shape  function,  as  it  flattens 
out  the  lateral  shape  function  near  beam  winds. 


Yaw  Moment  Due  to  Wind 


The  yaw  moment  coefficients  and  shape  function  are  also  strongly 
dependent  on  the  vessel  group  and  the  relative  longitudinal  location  of 
the  superstructure.  The  wind  yaw  moment  is  given  by  the  following 
equation: 

\  =  1/2  “a  \  L  CM  W  (A-M> 

where  =  wind  yaw  moment  coefficient  amplitude 

fM(ew)  =  normalized  yaw  moment  shape  function 

Reference  1  did  not  attempt  curve  fits  for  available  yaw  moment 
data.  However,  for  the  purposes  of  this  program  the  curves  that  appeared 
in  that  report  were  fitted  with  two  half  sine  waves.  Therefore  the 
shape  function  has  the  following  form: 

fM(0  )  =  sin  ♦  (A-25) 

n  w 

where  <t>  is  the  function  of  incident  wind  angle  and  vessel  group. 

The  variables  4>  and  CM  change  with  the  vessel  group  and  the  relative 
longitudinal  location  of  tne  superstructure.  The  method  used  for  deter¬ 
mining  wind  yaw  moment  is  not  as  analytical  as  the  wind  induced  longi¬ 
tudinal  and  lateral  loads.  However,  it  does  express  the  trends  found  in 
experimental  data. 
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Appendix  B 

BLOCK  FLOW  CHART,  SAMPLE  OUTPUT,  AND  LISTINGS 


STATMOOR  was  written  in  the  BASIC  language  for  Digital  Equipment 
Corporation's  (DEC)  Modular  Instrument  Computer  (MINC) .  This  version  of 
BASIC  is  typical  for  most  desk  top  computers.  If  the  program  is  to  be 
converted  for  use  on  another  computer,  the  following  notes  will  be  of 
value: 

1.  If  the  OR  and  AND  logical  operators  can  be  used  in  IF  statements  the 
program  can  be  made  significantly  shorter  in  some  sections. 

2.  The  MINC  BASIC  graphics  was  written  to  facilitate  data  acquisition. 
The  graphics  program  (MGRAPH)  will  need  to  be  rewritten  for  most  versions 
of  BASIC.  Shapes  were  defined  then  rotated  with  the  methods  explained 

in  Reference  13. 

3.  Because  of  the  limited  programming  work  space,  approximately 
7,000  words,  STATMOOR  was  broken  into  five  programs.  Except  for  slower 
interaction,  the  user  is  not  aware  of  program  segmentation.  The  COMMON 
and  CHAIN  statements  were  used  to  segment  the  program.  The  COMMON 
statement  appears  toward  the  beginning  of  each  program  and  protects  the 
listed  variables  when  each  program  is  run.  The  CHAIN  statement  replaces 
a  program  in  the  work  space  with  the  designated  program  and  continues 
execution  at  the  specified  line  number.  The  block  flow  chart  shown  in 
Figure  B-l  illustrates  the  relationship  between  the  programs. 

A  sample  tabular  output  is  shown  in  Figure  B-2,  which  summarizes 
the  user-supplied  information,  the  various  coefficients  and  final  magni¬ 
tudes  of  the  components  of  the  loads  at  equilibrium,  and  the  equilibrium 
conditions. 

A  listing  of  the  program  is  presented  as  Figure  B-3. 
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MOOR 

(main  program) 


Title  Page  — 


MEXPL 

(explanation  program) 


Master  Menu 

1.  Vessel  Input 

2.  Environmental  Input 

3.  Vessel  Review/Edit 

4.  Environmental  Review/Edit 

5.  Calculate - 

6.  Tabular  Output - 

7.  Graphic  Display - - 

8.  Exit  Program  —  j 


MCALC 

(calculation  program) 


MOUT 

(tabular  output  program) 


MGRAPH 

(graphic  display  program) 


Figure  R-l .  STATMOOR  block  flow  chart  showing  program  interactions, 


figure  B-2.  Sample  tabular  output  from  STATMOOR. 
PROGRAM  STATMOOR  -  STATIC  MOORING  ANALYSIS 
18-N0V-81  10149:34 

70  KDNT  TANKER  NSMB  TEST  *7561 


USER  INPUT 


VESSEL  INPUT 


VESSEL  GROUP!  TANKER 

Cl  TANKER  FULLY  BALLASTED  SUPERSTRUCTURE  LOCATION  CLUTTERED  DECK 

N  Y  STERN  N 


DISPLACEMENT  WATERLINE  LENGTH  BEAM  DRAFT 

87.600  LONG  TONS  740.00  FT  112.60  FT  43.90  FT 


WIND  AREA  INPUT  (AREAS  IN  SQUARE  FEET) 


END 


HULL  SIDE  SUPERSTRUCTURE  SIDE 


6.560  12.600 

MAX  HEIGHT  ABOVE  HULL 


6.420 

SUPERSTRUCTURE 


WATER  LINE  17.00  FT  79.00  FT 


ENVIRONMENTAL  INPUT 


CURRENT 

SPEED 

INCIDENT  ANGLE 

WATER  DEPTH 

4.00  KNOTS 

0.00  DEGREES 

200.  FT 

WAVE 

SIGNIFICANT  HEIGHT 

PERIOD  KNOWN 

SIGNIFICANT  PERIOD  INCIDENT  AN6LE 

12.00  FT 

Y 

7.40  SEC.  90.00  DEGREES 

WIND 

SPEED 

INCIDENT  ANGLE 

30.00  KNOTS  AT  33.0  FT  ABOVE  THE  WATER  90.00  DEGREES 


OUTPUT 

COEFFICIENT  NAXIMUMS 


CXB 

CXS 

THETAZ 

CY 

CMB 

CHS 

THETAM 

WIND 

0.700 

0.600 

100. 

0.725 

0.025 

0.110 

75. 

CURRENT  FRICTION 

0.002 

0.002 

FORM 

0.100 

0.100 

0. 

-.329* 

0.079* 

PROP 

1.000 

1.000 

*  -  INCLUDES  ANGLE  FUNCTION.  DEPTH/DRAFT  FACTOR 

=  1.051 

ENVIRONMENTAL  FORCES  RELATIVE  TO  VESSEL  AT  EQUILIBRIUM 


LOCAL 


FX  LBS 

FT  LBS 

M  FT-LBS 

ANGLE  (1 

CURRENT 

FRICTION 

6.240. 

FORM 

14.274. 

-505,678. 

89,662,400. 

PROP 

14,394. 

TOTAL 

34.907. 

-505,678. 

89,662,400. 

-30.45 

WAVES 

(MEAN) 

4,673. 

232,303. 

-2,761,260. 

59.55 

WIND 

8.341, 

40,270. 

-649,242. 

59.55 

STATIC  DESIGN  RESULTS 

SHIP  ANGLE 
30.4  DEGREES 
18 


HAWSER  LOAD 
237.979.  LBS 


HAUSER  ANGLE 
312.1  DEGREES 


MOOR,  main  program 


Figure  B-3.  STATMOOR  listing. 

to  m  mem  *statmoor*  static  mooring  analysis 

IS  REN  BY  JANES  V.  COX  CODE  LA A  CEL 

20  REN  26  AUGUST  1981 

21  ren  main  PROGRAM  ummmmtmitmmmmmtmm* 

22  REN 

23  REN 
2A  REN 

23  COMMON  F(19)>T<6)»C(ll>iH(5)>U(3!.B(3)»A(5)iV<A).H<3) 

30  COMMON  W(6.2).P(6.2).8(A.1).DA!7).RA(1) 

35  COMMON  RlfPl>T1.01»Ll»Sl.Hl.BltHO)EO»El»E2.NO.MliTO 
90  GO  TO  9300 
1600  CHAIN  'MCALC.BAC' 

5000  REM  VESSEL  CHARACTERISTICS  INPUT  STMT/ROUTINES  tmmMttmtttMt tttt 
5010  PRINT  ‘WHICH  OF  THE  FOLLOWING  VESSEL  GROUPS  DOES  THE  VESSEL  FIT  IN  ?* 

5020  PRINT  .‘1.  CARRIER* 

5030  PRINT  .‘2.  CRUISER* 

50A0  PRINT  .*3.  DESTROTER* 

5050  PRINT  .‘A.  CARGO* 

5060  PRINT  .‘5.  TANKER* 

5070  PRINT 

5080  GOSUB  5170  \  REM  GETS  CHOICE 

5083  IF  D10X1  GO  TO  5080 

5086  IF  D(0»5  GO  TO  5080 

5090  0(1)=0<0)  \  DA<1)=*N*  \  REM  NO  Cl 

5100  ON  D(l)  GOSUB  5120.5130.51A0.5150.5160 

5110  RETURN 

5120  M(0)=*CARRIER*  \  DA(3)=*S*  \  DA(5>=*N*  \  RETURN 

5130  DA(0)=*CRU1SER*  \  DA(3)=*D*  \  D*(5)=*M*  \  RETURN 

51  AO  DA(0)=*D£STROYER*  \  DA(3)-*D*  \  DA(5)=*N*  \  RETURN 
5150  DA(0)S*CARGO*  \  D(2)-3  \  D*<;A)=*M/A*  \  RETURN 
5160  DA(0)=*TANKER*  \  DA(3)=*S*  \  RETURN 

5170  PRINT  \  PRINT  \  PRINT  ‘ENTER  NUMBER  OF  ABOVE  CHOICE* i  \  INPUT  D'O)  \  RETURN 
5180  ERASE- TEXT ('TEXT')  \  MOVE-CURSOR T 1 f 1 )  \  RETURN 

5185  PRINT  ‘DESCRIBE  THE  VESSEL  (LIMIT  TO  15  CHARACTERS)’.  *)  \  INPUT  BA(O)  \  RETURN 
5190  PRINT  *IS  THE  VESSEL  A  CENTER  ISLAND  TANKER  (Y/N)*i  \  INPUT  DA(l)  \  RETURN 
5200  PRINT  *IS  THE  TANKER  FULLY  BALLASTED  (Y/N)*)  \  INPUT  DA(2)  \  RETURN 

5210  PRINT  ‘DOES  THE  VESSEL  HAVE  A  DISTRIBUTED  OR  SINGLE  SUPERSTRUCTURE  <D/S>*F  \  INPUT  D*<3>  \  RETURN 

5220  PRINT  ‘WHICH  OF  THE  FOLLOWING  BEST  DESCRIBES  THE  LOCATION  OF  THE  VESSEL'S  SUPERSTRUCTUR 

E  ?* 

5230  PRINT  ,*1.  STERN* 

52A0  PRINT  ,*2.  SUPERSTUCTURE  JUST  AFT  OF  MIDSHIP* 

5250  PRINT  ,*3.  CENTERED  SUPERSTRUCTURE- 

5260  PRINT  » *  A.  SUPERSTRUCTURE  FORWARD  OF  MIDSHIP* 

5263  PRINT  ,*5.  BOW* 

5270  PRINT 
5280  GOSUB  5170 

5283  IF  DT0K1  THEN  GO  TO  5280 
5286  IF  D(0)>5  THEN  GO  TO  5280 
5290  D(2)=D(0) 

5300  ON  D(2)  GOSUB  5310.5320.5330.53A0.5345 
5305  RETURN 

5310  DA(A)s*STERN*  \  RETURN 

5320  MCA) = ’SUPERSTRUCTURE  JUST  AFT  OF  MIDSHIP'  \  RETURN 
5330  DA( A) -'CENTERED  SUPERSTRUCTURE*  \  RETURN 
53A0  DA ( A ) = ' SUPERSTRUC  TURE  FORWARD  OF  MIDSHIP*  \  RETURN 
53A5  DA ( 4  > = *  BOW  *  \  RETURN 

5350  PRINT  ’DOES  THE  VESSEL  HAVE  CLUTTERED  DECKS  (I.£.  MASTS .  BOOMS.  PIPES.  AND  OThEF  SUBSTANTIAL  0 
BSTRUCTIONS)  <Y/N>*( 

5351  INPUT  DA(5)  \  RETURN 

5360  PRINT  'WHAT  IS  THE  DISPLACEMENT  OF  THE  'DA(O)’  IN  LONG  TONS')  \  INPUT  D1  \  RLI'JRn 
5370  PRINT  'WHAT  IS  THE  WATERLINE  LENGTH  IN  FEET*!  \  INPUT  LI  \  RETURN 
5380  PRINT  ’WHAT  IS  THE  BEAM  IN  FEET*.’  \  INPUT  B1  \  RETURN 
5390  PRINT  ‘WHAT  IS  THE  DRAFT  IN  FEET’)  \  INPUT  T1  \  RETURN 
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5400  PRINT  'HHAT  IS  THE  END  PROJECTED  AREA  OF  THE  VESSEL  ABOVE  THE  HATERLINE  IN  SQUARE 
BY  A  HEAD  HIND)'; 

5401  INPUT  A(5)  \  RETURN 

54'  PRINT  'HHAT  IS  THE  SIDE  PROJECTED  AREA  OF  THE  HULL  ABOVE  THE  HATERLIME  IN  SQUARE  FEET 
HIND)*?  \  INPUT  A(2> 

5411  RETURN 

5430  PRINT  *HHAT  IS  THE  AVERAGE  HEIGHT  OF  THE  TOP  OF  THE  HULL  ABOVE  THE  HATERLINE  IN 
T  H<2>  \  RETURN 

5440  PRINT  *HHAT  IS  THE  SIDE  PROJECTED  AREA  Of  THE  SUPERSTRUCTURE  ABOVE  THE  HATERLINE  IN 

AS  SEEN  BY  A  KAN  HIND)'! 

5441  INPUT  A<3)  \  RETURN 

5450  PRINT  'UHAT  IS  THE  AVERAGE  HEIGHT  OF  THE  TOP  OF  THE  SUPERSTRUCTURE  ABOVE  THE  HATERLINE 
NPUT  H<3)  \  RETURN 

5500  REN  ENVIRONNENTAL  CONDITIONS  INPUT  STMT/ROUTINES  ************************** 

5505  PRINT  'DESCRIBE  THE  MOORING  (LINIT  TO  15  characters):  •;/  input  b$U)  /  return 
5510  PRINT  , 'SPEED  IN  KNOTS  =  '»  \  INPUT  V(N)  \  RETURN  \  REN  CURRENT  OR  HIND 

5515  PRINT  , 'HATER  DEPTH  IH  FEET  =  'i  \  INPUT  HI  \  RETURN 

5520  PRINT  .'DIRECTION  IN  DEGREES  =  ')  \  INPUT  U(M)  \  RETURN  \  REN  ALL  LOADS. 

5530  PRINT  i 'SIGNIFICANT  HAVE  HEIGHT  IN  FEET  ^  ')  \  INPUT  HO  \  RETURN 

5540  PRINT  .'DO  YOU  KNOU  THE  SIGNIFICANT  HAVE  PERIOD  (Y/N)')  \  INPUT  D*(6)  \  RETURN 

5550  PRINT  .'SIGNIFICANT  HAVE  PERIOD  IN  SECONDS  =  '}  \  INPUT  TO  \  RETURN 

5560  PRINT  .'HEIGHT  OF  VELOCITY  NEASURENENT  IN  FEET  (E.G.  33  FT)  =  ')  \  INPUT  H(l!  \  RETURN 

5600  REN  VESSEL  CHARACTERISTICS  INPUT  CONTROL  ROUTINE  ttmmutmttmmtl 

5605  GOSUB  5180 

5607  PRINT  , 'VESSEL  CHARACTERISTICS  INPUT* 

5608  PRINT  Li  \  PRINT 

5609  GOSUB  5185  \  PRINT 

5610  GOSUB  5010  \  REN  VESSEL  GROUP 
5615  PRINT 

5620  ON  D(l)  GOSUB  5660.5660.5660,5700,5740 
5630  FOR  N=1  TO  9 
5635  PRINT 

5640  ON  N  GOSUB  5360,5370,5380,5390,5400,5410,5430,5440,5450 
5650  NEXT  N  \  GO  TO  9000 

5660  REN  WARSHIP  -  NO  ADDITIONAL  OBJECTIVE  INPUT  ttmmttmmmmtMt 
5670  RETURN 

5700  ren  cargo  ********************************** 

5710  GOSUB  5210 
5717  PRINT 

5720  IF  D$(3)-'S'  THEN  GOSUB  5220 
5725  PRINT 

5730  GOSUB  5350  \  RETURN 

5740  REN  TANKER  mmtttMtMtMttnmmttM 
5750  GOSUB  5190  \  PRINT  \  GOSUB  5200 
5755  B*(4)='N/A* 

5757  PRINT 

5760  IF  D*(1)-'N'  THEN  GOSUB  5220 
5765  PRINT 

5770  GOSUB  5350  \  RETURN 

5800  REN  ENVIRONNENTAL  CONDITIONS,  INPUT  CONTROL  ROUTINE  ******************** 

5810  GOSUB  5180 

5820  PRINT  , 'ENVIRONNENTAL  CONDITIONS  INPUT' 

5830  PRINT  L*  \  PRINT  \  PRINT 
5835  GOSUB  5505  \  PRINT 
5840  PRINT  'CURRENTS' 

5850  PRINT  ' _ ' 

5860  N=0  \  GOSUB  5510 
5865  GOSUB  5515 
5870  N=1  \  GOSUB  5520 
5880  PRINT 
5890  PRINT  'HAVES' 

5900  PRINT  ' _ ' 


FEET  (AS  SEEN 

(AS  SEEN  BT  A 

FEET')  \  INPU 
SQUARE  FEET  ( 

IN  FEET')  \  I 
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5910  GOSUB  5530  \  GOSUB  5540 
5920  T0=0 

5930  IF  D« ( 6 )  = ' Y *  THEN  GOSUB  5550 
5940  N=2  \  GOSUB  5520 
5950  PRINT 
5960  PRINT  ‘WIND* 

5970  PRINT  *_ _ * 

5980  GOSUB  5560 
5990  N-l  \  GOSUB  5510 
6000  N=3  \  GOSUB  5520 
6010  GO  TO  9000 

6100  REM  VESSEL  CHAR  REVIEW  STMT/ROUTINES  MtmmtmmttmmtMMt 
6105  PRINT  Nt'i  VESSEL  DESCRIPTION:  ,B»(0)  \  RETURN 
6110  PRINT  NV.  VESSEL  GROUP:  *BVO>  '  RETURN 
6120  PRINT  NV.  CENTER  ISLAND  TANKER:  *D$d:  \  RETURN 
6125  PRINT  NV.  FULLY  BALLASTED*.  ‘BSiO  '  RETURN 
6130  PRINT  NV.  SUPERSTRUCTURE  TYPE:  *  US  >  3  \  RETURN 
6140  PRINT  N4‘ .  SUPERSTRUCTURE  LOCATION:  *B* . 4 ■  RETURN 
6150  PRINT  NV.  CLUTTERED  DECKS:  *DV5>  \  REURh 
6160  PRINT  NV.  DISPLACEMENT  =  ,D1‘  LONG  TONS.1  \  f:ErURN 
6170  PRINT  NV .  WATERLINE  LENGTH  -  M*  FT.‘  \  RETURN 
6180  PRINT  NV.  BEAM  *  ’Bl*  FT.*  \  RETURN 
6190  PRINT  N$\  DRAFT  =  *T1*  FT.1  \  RETURN 
6200  PRINT  NV .  END  WIND  AREA  =  *A<5) ‘  SOUARf  FEET.*  \  RETURN 
6210  PRINT  NV.  HULL  SIDE  WIND  AREA  •*  ’A!2>*  SQUARE  FEET.*  \  RETURN 
6220  PRINT  NV .  TOP  OF  HULL  HEIGHT  ABOVE  WATERLINE  =  *H(2)*  FEET.’  \  RETURN 
6230  PRINT  Hi1.  SUPERSTRUCTURE  SIDE  WIND  AREA  =  *A<3>'  SQUARE  FEET.'  \  RETURN 
6240  PRINT  NV.  TOP  OF  SUPERSTRUCTURE  HEIGHT  ABOVE  WATERLINE  =  *H(3>*  FEET.’  \  RETURN 
6300  REM  ENVIRONMENTAL  CONDITION  REVIEW  STMT/ROUTINES  tittm 
6310  PRINT  *0.  RETURN  TO  MASTER  MENU.'  \  RETURN 
6315  PRINT  NV.  MOORING  DESCRIPTION.*  ’Bid)  \  RETURN 
6320  PRINT  NV.  CURRENT  SPEED  =  'V(O)'  KNOTS.*  \  RETURN 
6325  PRINT  N$*.  WATER  DEPTH  =  *H1*  FEET.*  \  RETURN 
6330  PRINT  NV,  CURRENT  DIRECTION  =  'U<1)*  DEGREES.  *  \  RETURN 
6340  PRINT  NV.  SIGNIFICANT  WAVE  HEIGHT  =  *H0*  FEET.*  \  RETURN 

6350  PRINT  N$‘ .  SIGNIFICANT  WAVE  PERIOD  KNOWN!  *W(6)  \  RETURN 

6360  PRINT  NV.  SIGNIFICANT  WAVE  PERIOD  =  'TO*  SECONDS.’  \  RETURN 

6370  PRINT  NV.  WAVE  DIRECTION  =  *U(2)*  DEGREES.*  \  RETURN 

6380  PRINT  N»\  WIND  VELOCITY  MEASURED  AT  *Hd)*  FEET  OFF  OF  WATER.*  \  RETURN 
6390  PRINT  NV ,  WIND  SPEED  =  *Vd)'  KNOTS.'  \  RETURN 
6400  PRINT  NV.  WIND  DIRECTION  =  *U<3>*  DEGREES,*  \  RETURN 
6480  B=ll  \  RETURN 
6485  B=14  \  RETURN 
6490  B=15  \  RETURN 

6500  REM  VESSEL  CHAR  REV/EDIT  CONTROL  ROUTINE 
6510  GOSUB  5180 

6520  PRINT  i* VESSEL  CHARACTERISTICS  REVIEW/EDIT* 

6530  PRINT  L*  \  PRINT  BVO>*  *Ptd)  \  PRINT  L»  \  PRINT 
6540  GOSUB  6310 

6550  ON  Dd)  GOSUB  6480,6480, *480)6485)6490 
6560  FOR  Nd  TO  B 
6570  NVSTRVN) 

6580  ON  Dd)  GOSUB  6700,6700,6700,6730,6760 
6590  NEXT  N 

6600  GOSUB  5170  \  REM  ASK  FOR  CHOICE 
6610  IF  D(0)"-0  THEN  GO  TO  90f") 

6620  IF  D\0»B  THEN  GO  TO  6600 
6630  ON  Dd)  GOSUB  6790.6790,6790,6820,4850 
6640  GO  TO  6500 

6700  rem  warship  review  mttmmmmtmtmmmt 

6710  ON  N  GOSUB  6105,6110, 6160, 6170, 6180, 6190,6200-6210,6220, f23<*  6240 
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6720  RETURN 

6730  rem  cargo  revieu  ttmmmmmmmmmm 

6740  ON  N  60SUB  6105,6110,6130,6140,6150,6160,6170,6180,6190,6200,6210,6220,6230,6240 
6730  RETURN 

6760  REM  TANKER  REVIEW  MUmMMttttttMtmtMt 

6770  ON  N  GOSUB  6105,6110,6120,6125,6140,6150,6160,6170,6180,6190,6200,6210,6220,6230,6240 
6780  RETURN 

6790  reh  warship  edit  itMtmmmtiUMttnmttnt 

6800  ON  0(0)  GOSUB  5185.50107 5360 * 5370 » 5380 f 5390 » 5400 f 5410» 5430 f 5440 . 5450 
6810  RETURN 

6820  REN  CARGO  EDIT  mmttMtnnMtmnttUttM 

6830  ON  D(0)  GOSUB  5185 » 5010.5210.5220.5350. 5360.5370. 5380. 5390 7  5400.5-410. 5430.5440.5450 
6840  RETURN 

6850  ren  tanker  edit  mtnmtmtMtMtmmutnt 

o860  ON  D(0)  GOSUB  5185.5010.5190.5200.522075350.5360.5370.5380.5390.5400.5410.5430.5440.5450 
6870  RETURN 

6900  REN  ENVIRONMENTAL  CONDITION  REVIEW/EDIT  CONTROL  ROUTINE  mtitmtUttUtl 
6910  GOSUB  5180 

6920  PRINT  ENVIRONMENTAL  CONDITION  REVIEW/EDIT 1 

6930  PRINT  L*  \  PRINT  Bi(0)*  *84(1)  \  PRINT  Li  \  PRINT 

6940  GOSUB  6310 

6950  FOR  N=1  TO  11 

6960  Ni=STR4(N) 

6970  ON  N  GOSUB  6315. 6320. 6325, 6330. 6340. 635C.e340i6370.6380.6390. 6400 
6980  NEXT  N 

6990  GOSUB  5170  \  REH  ASK  FOR  CHOICE 
7000  IF  D(0)<>0  THEN  GO  TO  9000 
7010  IF  D(0)>11  THEN  GO  TO  6990 

7020  ON  D(0)  GOSUB  5505,7040,5515,7050,5530,5540,5550,7060,5560,7080,7090 

7030  GO  TO  6900 

7040  N=0  \  GOSUB  5510  \  RETURN 

7050  N=1  \  GOSUB  5520  \  RETURN 

7060  N=2  \  GOSUB  5520  \  RETURN 

7080  N=1  \  GOSUB  5510  \  RETURN 

7090  N=3  \  GOSUB  5520  \  RETURN 

7230  PRINT  \  PRINT  \  PRINT  \  PRINT  ,, ‘PRESS  RETURN  TO  CONTINUE.1!  \  INPUT  Bi(7)  \  RETURN 
8200  CHAIN  'NOUT.BAC' 

9000  ren  ********  master  menu  unttm 
9010  ren  ******************************** 

9015  Li** _ 

9020  GOSUB  5180 
9030  PRINT  , . ’MASTER  MENU’ 

9040  PRINT  Li  \  PRINT  \  PRINT 
9050  PRINT  ,’l.  INPUT  VESSEL  CHARACTERISTICS.* 

9040  PRINT 

9070  PRINT  ,‘2.  INPUT  ENVIRONMENTAL  CONDITIONS.1 
9080  PRINT 

9090  PRINT  ,‘3.  REVIEW/EDIT  VESSEL  CHARACTERISTICS,* 

9100  PRINT 

9110  PRINT  ,*4.  REVILii/EDIT  ENVIRONMENTAL  CONDITIONS.* 

9120  PRINT 

9130  PRINT  ,*5.  CALCULATE. * 

9140  PRINT 

9150  PRINT  ,’6.  TABULAR  OUTPUT,* 

9160  PRINT 

9170  PRINT  »*7,  GRAPHIC  DISPLAT.* 

9180  PRINT 

9183  PRINT  ,*8.  EXI’  PROGRAM.* 

9186  PRINT 
9190  GOSUB  5170 
9200  IF  0(0X0  GO  TC  9000 
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*nu  ip  utoi>»  eo  to  woo 

9220  OX  D(0)  GO  TO  5600 r 5800 » 0500 , 6900 »1 600 , 8200 il 0500 t 9230 
9230  GQSU8  5180 
9240  PRINT  'BYE  BYE* 

9250  STOP 

9300  REH  TITLE  SCREEN  mMttUtmtMM 

9305  A$=  * ' CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC  * 
9310  60SUB  5180 

9320  N0VE.CURS0R<3»1)  \  GOSUB  9495  \  GOSUB  9490 
9335  PRINT  USING  A4»**  STATHOOR  ** 

9340  GOSUB  9490  \  GOSUB  9495  \  PRINT 
9355  PRINT  USING  A4» 'STATIC  HOORING  ANALYSIS* 

9360  PRINT  USING  A4,*V£RS10N  1.0* 

9365  PRINT  using  a$/april  198T 
9370  MOVE-CURSOR < iSf i ) 

9375  PRINT  USING  A4»‘BY  JANES  0.  COX* 

9385  PRINT  USING  A*. 'CODE  L44* 

9387  PRINT  USING  A4»*NAVY  CIVIL  ENGINEERING  LAB* 

9390  PRINT  USING  A4»*P0RT  HUENEHEt  CALIFORNIA  93043* 

9415  FOR  N=1  TO  3000  \  NEXT  N 
9420  GO  TO  9500 

9490  PRINT  USING  A4.*l  *'  \  RETURN 

9495  PRINT  USING  Ati'limmitmt*  \  RETURN 
9500  REH  EXPLANATION/ ASSUNPTION  PAGE  tmtttt 
9510  GOSUB  5180 

9520  PRINT  *D0  YOU  WISH  TO  REVIEW  THE  PROGRAM  EXPLANATION  (Y/H)*i  \  INPUT  DS(7> 

9530  IF  M(7)0*N*  T1£N  CHAIN  'HEXPL.BAC' 

9550  GO  TO  9000 
10500  CHAIN  'HGRAPH.BAC' 

32000  END 
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10  REM  PROGRAM  •STATHOOft*  STATIC  H00RING  ANALYSIS 
15  REM  BY  JANES  V.  COX  CODE  L44  CEL 

20  REM  26  AUGUST  1981 

21  REN  EXPLANATION  PROGRAN  tUtWmtmUMtUMttmmttUmtU 

22  REN 

23  REN 

24  REN 

25  CONNON  F(19)»T(6).C(U)»N(5)»U<3).D(3)tA(5).V(4).H<3) 

30  CONNON  M(6,2),P(6»2).B(4fl)»M(7)»B$(l) 

35  CONNON  RliPl.Tl>Hl>Ll,Sl.Hl»Bl,H0.E0,El,E2,K0»Ml.T0 
50  GOSUB  5180 

60  PRINT  ’  PROGRAN  STATMOOR  -  STATIC  NOORING  ANALYSIS’ 

70  PRINT  \  PRINr  ’OVERVIEW*  \  PRINT  * - ’ 

80  PRINT  ’THIS  PROGRAN  CALCULATES  THE  EQUILIBRIUM  CONDITION  OF  A  SINGLE  POINT  NOOSED’ 

90  PRINT  ’SHIP  ACTED  UPON  BY  STEADY  CURRENT.  WAVE,  AND  WIND  LOADS.  THE  LOADING  CURVES’ 

100  PRINT  'HERE  OBTAINED  FROM  A  VARIETY  Of  SOURCES.  AND  NO  DYNAMICS  HERE  INCLUDED.’ 

110  PRINT  ’FOR  FURTHER  INFORMATION  SEE  REPORT  TN44-82-1.’ 

120  PRINT  \  PRINT  ’PROGRAN  OPERATION’  \  PRINT  • - ’ 

130  PRINT  ’THE  PROGRAN  IS  WRITTEN  IN  A  CONVERSATIONAL  NODE.  YOU  NUST  RESPOND  T0‘ 

140  PRINT  ’QUESTIONS  AND  CHOICES.  INPUT  CONSIST  OF  TYPED  DESCRIPTIONS.  NUMERICAL* 

150  PRINT  ’VALUES.  AND  LETTERS.’ 

160  PRINT  ’*  TYPED  DESCRIPTIONS  -  AT  ONE  POINT  IN  THE  PROGRAN  YOU  HILL  BE  ASKED  TO' 

170  PRINT  ’DESCRIBE  THE  SHIP.  FOR  EXANPLE  YOU  HA>  TYPE  70  KLWT  TANKER  FOLLOWED* 

180  PRINT  'BY  THE  RETURN  KEY.’ 

190  PRINT  ’*  NUNERICAL  VALUES  -  THEY  NUST  BE  TYPED  IN  HITH  NO  ALPHA  CHARACTERS  AND’ 

200  PRINT  ’FOLLOWED  BY  THE  RETURN  KEY.’ 

210  PRINT  ’*  LETTERS  -  AT  SEVERAL  POINTS  THE  PROGRAN  ASK  A  QUESTION  CONTAINING  THE’ 

220  PRINT  ’CHARACTERS  (Y/N).  1  HIS  MEANS  YOU  SHOULD  ENTER  Y  FOR  YES  OR  N  FOR  NO,’ 

230  PRINT  ’FOLLOWED  BY  THE  RETURN  KEY.’ 

240  PRINT  \  GOSUB  7230 
250  GOSUB  5180 

260  PRINT  ’PROGRAM  OPTIONS’  \  PRINT  ’ . - . ’ 

270  PRINT  ’»  INPUT  SECTIONS  -  YOU  CAN  ENTER  ACL  THE  INFORMATION  ABOUT  THE  VESSEL* 

280  PRINT  ’OR  THE  ENVIRONMENTAL  LOADS-' 

290  PRINT  ’»  REVIEW/EDIT  SECTIONS  -  YOU  CAN  REVIEW  AND  EDIT  INDIVIDUAL  VALUES  OF  THE’ 

300  PRINT  ’INPUT  FOR  ThE  VESSEL  OR  THE  ENVIRONnENAL  LOADS.’ 

310  PRINT  ’*  CALCULATION  SECTION  -  CALCULATES  EQUILIBRIUM  CONDITION  ;WITH  O  OPTION  TO’ 

315  PRINT  ’FIND  THE  WORST  CURRENT  ANGLE)  AND  HAS  A  BRIEF  VIDEO  OR  PRINTER  OUTPUT • 1 
320  PRINT  ’»  TABULAR  OUTPUT  SECTION  -  YOU  CAN  OBTAIN  VIDEO  SR  PRINTER  OUTPUT.’ 

330  PRINT  ’*  GRAPHIC  OUTPUT  SECTION  -  PLOTS  EQULIBRIUN  CONDITIO-  OF  SMIP  AnD  HRlSE1’ 

340  PRINT  ’WITH  INCIDENT  ANGLES  OF  ENVIRONMENTAL  LOADS.’ 

350  PRINT  ’t  EXPLANATION  -  YOU  CAN  USE  ThE  GRAPHICS  DEMO  TO  CREATE  GRAPHIC  OUTPUT ’ 

360  PRINT  ’INDEPENDENT  CF  THE  CALCULATED  RESULTS.’ 

370  PRINT  \  PRINT  ’ANGLE  CONVENTION’  \  PRINT  . . ’ 

380  PRINT  ’THERE  ARE  TWO  ANGLE  CONVENTIONS  IN  THE  PROGRAM,  ABSOLUTE  A'<D  LOCAL  • 

3°0  PRINT  ’t  ABSOLUTS  -  THIS  ANGLE  IS  ZERO  AT  TWELVE  O'CLOCK  On  THE  VIDEO  an;  INCREASES’ 

400  PRINT  ’COUNTERCLOCKWISE.  IT  IS  USED  TO  ENTER  TH£  E“'IRGNMENTAL  LOAD  PIPESVM’ 

410  PRINT  ’AND  OUTPUT  THE  SHIR  AND  HAWSER  ANGLES.' 

420  PRINT  •«  LOCAL  -  THIS  ANGLE  IS  RELATIVE  TO  THE  VESSEL  DIRECTION.  HEAD  CASING’ 

430  PRINT  ’IS  ZERO  DEGREES  AND  THE  ANGLE  INCREASES  COUNTERCLOCKWISE.  IT  IS  UiSD  r0F 1 

440  PRINT  'CALCULATION  OF  LOADS  AND  OUTPUT  IN  THE  FORCE  TA8L£-’ 

450  GOSUB  7230 
460  GOSUB  5180 

470  PRINT  ’LOCAL  COORDINATE  SYSTEM’  \  PRINT  ’ - ’ 

480  PRINT  ’ALL  ENVIRONMENTAL  LOADS  IN  THE  OUTPUT  TABLE  ARE  IN  THE  LOCAL  COO‘D:',p!S  S*CTEN.’ 
490  PRINT  ’*  X  -  rOSITIVE  AT  THE  LOCAL  ANGLE  OF  180  DEGREES.  TOWARD  THE  STERN. ‘ 

500  PRINT  ’*  V  -  POSITIVE  AT  THE  LOCAL  ANGLE  OF  270  DEGREES,  TOWARD  THE  STARBOARD. ' 

510  PRINT  ’*  M  -  POSITIVE  COUNTERCLOCKWISE.’ 

520  PRINT 

530  PRINT  ’GRAPHICS  DEMO'  \  PRINT  ’ - ’ 

540  PRINT  ’THE  FOLLOWING  DEMONSTRATION  IS  TO  FAMILIARIZE  ThE  USER  WITH  THE  ABSOLUTE’ 
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550  PRINT  'ANGLE  CONVENTION,  THE  ANGLE  ENTERED  FOR  THE  ENVIRONMENTAL  LOADS  IS  THE' 
560  PRINT  'INCIDENT  ANGLE,  THESE  LOADS  ARE  REPRESENTED  AS  FOLLOWS'.' 

565  PRINT  *C  *  CURRENT'  \  PRINT  'D  =  NAVE  DRIFT  FORCE'  \  PRINT  '«  *  KIND'  \  PRINT 

566  PRINT  \  PRINT  'THE  FIRST  EXAMPLE  IS  FOR  ALL  ANGLES  AT  ZERO  DEGREES.' 

567  PRINT  'ADDITIONAL  EXAMPLES  WITH  USER  DEFINED  ANGLES  FOLLOW.'  \  PRINT 
570  GQSUB  7230 

575  DISPLAY-CLEAR 

580  B»(0)='GRAPHICS*  \  M(1)='DEM0*  \  Dtt7)='DEH0' 

585  H0=1  \  V(0>"1  \  V(1M 
650  CHAIN  'M6RAPH.BAC' 

5180  ERASE-TEXT ('TEXT')  \  MOVE-CURSOR (1,1)  \  RETURN 

7230  PRINT  # ‘PRESS  RETURN  TO  CONTINUE'?  \  INPUT  D*<7>  \  RETURN 

32000  END 


MCALC,  calculation  program 


10  REM  PROGRAM  •STATHQOR*  STATIC  MOORING  ANALYSIS 
IS  REN  BY  JANES  V.  COX  CODE  L44  CEL 

20  REN  26  AUGUST  1981 

21  REN  CALCULATION  PROGRAM  tttUttttitttttt 

22  REN 

23  REN 

24  REN 

25  COMMON  F(19)»T(6),C(11)>H(5)»U(3)»D(3j,Ai5),V(4),H(3) 

30  COMMON  y<6,2),P(6>2),B(4<l)»D4(7),B4(l) 

35  COMMON  R1,P1,T1,D1»L1»S1»H1,B1»HO>£0>£1>E2»HO,N1,TO 
40  DEF  FND(X)*X»180/PI  \  REN  CONVERSION  FROM  RAD  TO  DES. 

43  DEF  FNR(X)=XtPI/180  \  REN  CONVERSION  FROM  DES  TO  RAD 

45  ren  fourier  cocf.  DATAmmmmmm 

50  DATA  .908,0,-. 116,0, -.033 

53  DATA  -.0252.-.0904..0032.. 0109,. 0011 

55  REN  DRIFT  FORCE  TRANSFORM  DATA 

57  REM  FX  LAMDA/l  TRANS 

60  DATA  1,. 65,. 4,. 25 

63  DATA  0,0,. 07,. 08,0 

65  REN  FY  KT  TRANS 

67  LATA  .25,. 75, .97,3 

70  DATA  0, .01, .64, .81,1 

73  REN  N  LANDA/L  TRANS 

75  DATA  1.4, .95, .6, .25 

77  DATA  0,0, -.02, -.02,0 

78  03=0  \  GOSUB  5180  V  PRINT  'DO  YOU  UISH  TO  FIND  THE  CURRENT  DIRECTION  (WITHIN  10  DEGREES'  WHICH  PRODUCES* 

79  PRINT  ‘MAXIMUM  HAUSER  LGAD  <Y/N>  *S  \  INPUT  014 

80  PRINT  ‘DO  YOU  WANT  VIDEO  OR  PRINTED  RESULTS  V/P)  •?  \  INPUT  024 

81  02=0 

32  IF  024= *P*  THEN  02=1 

83  OPEN  'LP1'  FOR  OUTPUT  AS  FILE  1 

84  GO  TO  1600 

95  PRINT  102,,, 'CALCULATION  OUTPUT* 

96  L4=* _ _ * 

87  PRINT  102, L4  \  PRINT  #02,64(0) *  *B4!1)  \  PRINT  402, Li  \  PRINT  402, 

88  PRINT  402, 'HAUSER  LOAD  HAWSER  ANGLE  SHIP  ANGLE* 

89  A4=*4444,444  LBS  4444.4  DEGREES  4444-4  DEGREES*  \  RETURN 

95  PRINT  402, USING  A4,F(0),T(6),E2  \  RETURN 

100  REN  FX  CURRENT  SUBR0UTINE«mm*4 

110  REN  CONSIDER  FORM, FRICTION, AND  PROPELLER  DRAG 

120  REN  REF.  DN-26  OTC  ALTMAN 

125  IF  V(0)=0  THEN  GO  TO  110 

130  A=COS(B)  \  R2=ABS(R1*A)  \  C <1 > = . 075/ ( LOSlO < R2  — 2  *2 

140  F(2)=F(t)*A  \  REN  PROP 

150  F(4)=F(3)*A"3*C(1)  \  REM  FRIC 

160  F(6)=F(5)4A“34.1  \  REN  FORM 

170  F(7)=F(2HF(4HF<6> 

180  RETURN 

200  REN  FOURIER  SERIES  COEF  CALC '£»♦** 

210  C=0 

220  FOR  14=0  TO  4  \  A=N+1  \  C=C+B'N,D>»S!N'A*Bi  \  NEXT  N 
230  RETURN 

250  REN  FY  CURRENT  SUBR0UT!ME»»mm*»4m 
260  D=0  \  SOSUB  200  \  REN  DETERMINES  CT-B> 

270  C(3i=C  \  F(8)*F(9itC(3; 

290  RETURN 

350  REN  N  CURRENT  SUBROUT  INE44WM 
360  D*1  \  GOSUB  200  \  REN  DETERMINES  th(B; 

370  C(11)=C  \  H(0»=M(1)4C'.11.> 

380  RETURN 

400  REM  FX  DRIFT  SUBROUTINE#**#* 
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410  F<11)=F(10>*CGS<B>*ABS<CGS<BM 
420  RETURN 

450  REN  FT  DRIFT  SUBROUTINE****! 

440  F(13)=F(i2)*SlNiB)ttBS(SlN;B>» 

470  RETURN 

500  REN  «  DRIFT  SUBROUTINE************ 

510  H<3)=N(2)*S1N(2*B) 

520  RETURN 

650  REN  DEFINE  THE  EXCITATION  AT  CONSTANT  U  INTERVALS  AND 
455  REN  INTERPOLATE  THE  TRANSFER  FUNCTION  TO  GET  THE 
660  REN  RESPONSE 

670  REN  USE  F  AS  THE  INSTANTANEOUS  FREQUENCE  VALUE. 

680  U=10 

690  IF  H0>25  THEN  Lf=  14 

700  F=PI/‘J 

710  B=0 

720  K=2 

730  D=PI/60 

740  FOR  1=1  TO  60 

759  2=01/F"4 

760  IF  Z>25  THEN  E=1.00000E-}6 
770  IF  Z>25  GO  TO  790 

780  E=A1/F"5*EXP(-Z) 

790  IF  F<=«U..N»  30  TO  S20 
800  K=K+1 
810  GO  TO  790 

820  ft=P<K-l»N»+(F-U(K-l »N>  > '1U(K»N)-U<K-1»H> )*(P(K»N)-P(N-1tN! > 
830  R=E*A 
840  B=B+R*D 

850  IF  FC2.6  GO  TO  990 

860  IF  D=PI/60  THEN  f=F+D/2-t*2 

870  REN  NAKE  NEW  FREQUENCY  INTERVAL  *  4*011  INTERVAL 

880  IF  D-PI/60  THEN  D=D<4 

890  F=F+D 

900  NEXT  I 

910  RETURN 

950  REN  FX  WIND  SUBROUTINE  ****** 

960  REN  DETERMINE  SHAPE  FUNCrION 

970  IF  D(3)=0  THEN  GO  TO  1050 

980  REN  DEFAULT  DISTRIBUTED  SS 

990  IF  ABS<K)  :=T THEN  Z«A8i(!'-J»90/T*.4-490 

1000  IF  ABSUvttT',4)  THEN  Z*ABS,-F  *(?•)■' f  18C-T (i;i )tl90-90*T '(’3 

1010  REN  CONVERT  TO  RADIANS 

1020  Z=FNR',Z; 

1030  A=(SIN<Z'-SIN'5*Z)/10) '.9 

1040  SO  TO  1100 

1050  REN  SINGLE  DISTINCT  SS 

1060  IF  ABS(KK=T(*>  THEN  ;=APS(K'*90/T<*> 

1070  IF  ABS (K )>T ( 4 >  TNEN  >ABg!*>-M;  »90  ••  18C-S .  »>)+«•) 

1090  Z=FNR(Z> 

1090  A=COS(Z) 

1100  REN  DETERNINE  CQ£r.  CXB  OP  CXS 
1110  IF  AB8<KK=T<4.-  THEN  C(7;=C<4> 

1120  IF  ABS(K»T(4)  THEN  C‘.7:=C't; 

1130  F(15)=F'-14)*A*C(7) 

1140  RETURN 

1200  REN  FY  WIND  SUBROUTINE********** 

1210  Z-FNRtK; 

1220  A=(SIN-,Z)-SINi5*Z?/20'/  1-1/20) 

1230  F<17:=C,;:6;*A 
1240  RETURN 
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1300  REM  H  KIND  SUBRQUTINEammW 
1310  ROI  DETERMINE  SHAPE  FUNCTION 
1320  IF  ABS(KX=T (5)  THEN  GO  TO  1^70 
1330  REN  THETA) THETA  2 
1335  Z=K-T(5)  \  IF  K<0  THEN  Z=K+T(5> 

1340  Z=Z*180/(180-T(5» 

1350  C(10)*C<9> 

1340  GO  TO  1400 
1370  REM  TH£TA<=TH£TA  2 
1380  Z*K*180/T(5) 

1390  C(10)=-C(8) 

1400  REM  CONVERT  TO  RADIANS 
1410  Z=FNR<Z> 

1420  H(5)*M(4)*C(10)*SIN(Z1 
1430  RETURN 

1435  FOR  05=0  TO  350  STEP  10  \  d<l»=C5  \  GO  TO  1420 

1440  IF  U(1)=0  THEN  60SU6  5180 

1441  IF  U(1)*0  THEN  GOSUB  85 

1442  PRINT  W2f 'CURRENT  ANGLE  *  'U<1> 

1444  GOSUB  95  \  PRINT  402* 

1444  IF  F(0)<03  THEN  GO  TO  1450 
1448  03*Fl0)  \  04=U:U 

1450  NEXT  05 

1452  014=' LAST*  \  UU)*04  \  GO  TO  1420 

1440  GOSUB  85  \  GOSUB  95  \  GOSUB  7230  \  DI«7>=*N‘  \  GO  TO  1934 

1500  REM  CALCULATION  OF  SUMH)*ttM*m 

1503  IF  El)340  THEN  £l=£l-3*0 

1504  IF  E0>=340  THEN  E0=E0*340 

1505  FOR  N=l  TO  3  \  7<H)=U<N)-E1  \  NEXT  N 
1510  B=FNR(T(1)) 

1515  GOSUB  250  \  REfi  FT  C 
1520  GOSUB  350  \  REM  M  C 
1525  B=FNR<T(2)> 

1530  GOSUB  450  \  REM  FT  D 
1535  GOSUB  500  \  REM  M  D 
1540  K*r<3# 

1545  IF  T(3)>180  THEN  K=H-340 
1547  IF  K3X-1S0  THEN  K=K+3eO 
1550  GOSUB  1200  \  REM  FT  U 
1555  6QSUB  1300  \  REM  M  U 
1540  REM  SUM  MOMENTS 

1545  Ml=M(0)+M(3)+M(5)4Ll/2*'.F(8)tF(13)+F<17)) 

1570  RETURN 

1580  RBI  SECTION  OF  INTERP.  LOOP  IF  SOLUTION  IS  HIT  EXACTLY 

1585  Zl-4  \  GO  TO  1850 

1400  REM  CALCULATION  SECTIONWmWttmn 

1410  REM  PRELOOP  SECTION 

1415  60SUB  2000  \  REM  CALC  INIT  SUBROUTINE 

1417  IF  Ol^'Y'  THEN  GO  TO  1435 

1420  11*30  \  REN  ANGLE  INCREMENT 

1430  E1*0 

1440  GOSUB  1500  \  REM  ENTER  HERE  FOR  SECOND  SOLUTION. 

1445  REM  LOOP  SECTION 
1450  E0-E1 

1440  M0*B1 
1470  El*Eim 
1480  GOSUB  1500 
1490  D=NC*M1 

1700  IF  D>0  THEN  00  TO  1450  \  REM  CONTINUE  INCREMENT 
1710  IF  D<0  THEN  60  TO  1790  \  REM  STABLE  OR  UNSTABLE  SOLUTION 

1720  REN  IH) 

1730  IF  N0*0  GO  TO  1750 
1740  IF  N1«0  GO  TO  1740 
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1750  IF  HKO  THEN  SO  TO  1770 

1755  SO  TO  1650  \  REH  UNSTABLE 

1760  IF  H0>0  THEN  SO  TO  1780 

1765  60  TO  1650  \  REN  UNSTABLE 

1770  E2=£0  \  SO  TO  1860  \  REM  SKIP  ZERO  PROCESS 

1780  E2*E1  \  GC  TO  1860  \  REH  SKIP  ZERO  PROCESS 

1790  IF  MIXMl  GO  TO  1650  \  REN  UNSTABLE 

1800  REH  OTHERWISE  WE  HAVE  SOLUTION  ENVELOPE 

1805  REM  CONVERGE  ON  SOLUTION  WITHIN  ENVELOPE 

1810  REH  SET  ACCURACY  TO  BETTER  THAN  1  DEGREE  >  A  LOOPS 

1815  E2*C1  \  H2=H1 

1820  FOR  Zl=l  TO  4  \  El =E0+ ( E2-E0 ) *H0/ ( M0-M2 ) 

1830  60SUB  1500 

1835  IF  M1«0  THEN  60  TO  1580  \  REM  SOLUTION  AT  £1 

1840  IF  M1>0  THEN  £0=E1 

1842  IF  H1>0  THEN  H0=H1 

1845  IF  HKO  THEN  E2=E1 

1847  IF  HKO  THEN  H2*M1 

1850  NEXT  Z1 

1855  E2*£l  \  REH  SHIP  AN6L£=E2 

1860  FOR  N=1  TO  3  \  T(N)=U(N)-E2 

1861  IF  T(NKO  THEN  T(N)=T(N)+360 

1862  NEXT  N 
1865  B=FNR(T(1)1 

187©  GOSUB  100  \  REH  FX  C 
1875  GQSJ©  .750  \  REN  FY  C 
1878  GOSt®  350  \  REH  h  C 
1880  B=FNR(T(2>) 

1885  GOSUB  400  \  REN  FX  B 
1890  GOSUB  450  \  REH  FY  0 
1893  GOSUB  500  \  REH  M  0 
1895  K=T(3) 

1900  IF  T(3)>180  THEN  K*Tv3;-360 

1905  GOSUB  950  \  REH  FX  W 

1910  GOSUB  1200  \  REN  FY  W 

1913  GOSUB  1300  \  REH  H  W 

1915  F<18)=-CF<7)+F( U)+F<15»)  \  REH  FX  HAWSER 

1920  F(19)»-(F(8)+F(13)+F<17))  \  REH  FY  HAWSER 

1923  REH  DETERMINE  HAWSER  ANGLE 

1925  IF  F<18)=0  THEN  GO  TO  1990 

1930  T<6)=ATN(F(19)/F(18)) 

1940  IF  F(18»0  THEN  T(6t*T(6>m 
1950  T(6)=FND(T(6))4E2 
1955  IF  T(6)>«360  THEN  T(6)=Tl6!-360 
1960  IF  T(6K0  THEN  T(6)*n6)+360 
1970  F(0)«S8R(FU8r2+F(195',2) 

1972  IF  Olls'Y'  THEN  60  TO  1440 

1973  IF  01$**LAST'  THEN  SO  TO  1460 

1974  IF  D4(7)»*r*  THEN  60  TO  1980 

1975  GOSUB  5180 

1980  GOSUB  85  \  GOSUB  95 

1983  PRINT  \  PRINT  \  PRINT  *B0  YOU  WISH  TO  LOOK  FOR  ANOTHER  STABLE  SOLUTION  <Y/N)S  \  INPUT  D»<7 

1984  IF  WmO'Y*  THEN  CHAIN  'HOOR.BAC'  LINE  9000 

1985  El*( INT(E2/30)B1 )I30  \  IF  El*360  THEN  E1«0 

1986  60  TO  1640 

1990  IF  F(19)>*0  THEN  MM.5IPI 
1993  IF  F(19)<©  *HEN  T<6>«.5tPI 
1995  80  TO  1960 

2000  REH  CALCULATION  INITIALIZATION 

2011  BEX  SECTION  1...  BASED  ON  SHIP  TYPE  SETS  PROPERTIES  . 

2015  REH  E.6.  PROP  AREA.  MOMENT  CURVE  ETC. 


2020  REM  SET  FLAB  FOR  X  -  CURRENT  SHAPE  FUNCTION 
2030  D(3)=l 

2040  IF  BMD-'S'  THEN  D(3)=0 

2050  RED  ASSIGN  HIND  MOMENT  AND  FX  CURVE  VALUES. 

2060  K=D(1) 

2070  IF  K=5  T1CN  60SUB  2480 

2080  ON  K  GOSUB  2300. 2340 .2370. 2400 » 2425 i 2425 

2090  REM  PROPELLER  AREA  CALCStllimmitt 

2100  ON  DID  GOSUB  2240.2250.2260.2270.2280 

2110  A(0)=L1*B1/N  \  Pl=l 

2120  A(1>*A(0)/U.062-.229«P1) 

2130  REM  SPECTRUM  CHOICE 
2140  IF  M(6)='Y*  THEN  GOSUB  2220 
2150  IF  D4(6)**N‘  THEN  GOSUB  2230 
2160  GO  TO  2500 

2220  REM  BRETSCHNEIDER  HIND  UAVE  MODEL 
2222  IF  TO^O  THEN  GO  TO  222V 
2225  A1*262.5IH0“2/T0',4  \  01=1050/ T0~ 4 

2228  RETURN 

2229  Ql=9.90000E+33  \  RETURN 

2230  REM  PIERSON  MOSKOWITZ  HIND  HAVE  MODE: 

2233  IF  H0=0  THEN  GO  TO  2237 

2235  Al=8.4  \  Ql=33.56/H0"2  \  RETURN 
2237  Al=8,4  \  01=9.90000EF32  s  RETURN 
2240  N-125  \  REM  CARRIER  AREA  RATIO  TAR! 

2243  IF  LUBIOIOOOOO  THEN  N=28« 

2246  RETURN 

2250  N=160  \  RETURN  \  REM  CRUISER  AR 
2260  N=100  \  REM  DESTROYER  AR 
2263  IF  L1*B1<=12000  THEN  N=210 
2266  RETURN 

2270  N=240  \  RETURN  \  REM  CARGO  AR 

2280  (4=270  \  RETURN  V  REM  TANKER  AR 

2290  IF  M(5)*'N*  THEN  RETURN  \  REM  CLUTTERED  DECK  ADDER 

2295  C<4)=C«4H.08  \  C<6T=C!6K.08  \  RETURN 

2300  T ( 4)=120  \  C(4)=.4  \  C'6>=.4  \  REM  FX  CARRIER 

2310  T(5)=90  \  C(8)- .071  \  Cv9)=,Q71  \  REH  ft 

2315  IF  D(l)«5  THEN  C(4>=,7  \  REM  BAl  TANKER  CHANGES 

2317  IF  D(l)=5  THEN  C(6>%6 

2320  RETURN 

2330  T(4)=li0  \  C(4)*.7  \  C(6)*.8  \  RETURN  \  REM  FX 
2340  GOSUB  2330  \  REM  FX  CRUISER 
2350  T(5) =90  \  C(8)=.056  \  C<9)=.0o3  \  REM  M 
2360  RETURN 

2370  GOSUB  2330  \  REM  FX  DESTROYER 

2380  T (5>=110  \  C!8>=.123  \  C(9)=,018  \  REM  M 

2390  RETURN 

2400  F=D(2)  \  REM  CARGO 

2410  IF  M(3)=’D'  THEN  F=3  \  REM  DISTRIBUTED  -  CENTERED 
2415  GCSUB  2455 

2420  C( 4)= .7  \  C(6)s.6  \  GOSUB  2290  \  RETURN 
2425  F=D<2)  \  GOSUB  2455  \  REM  TANKER 
2430  C(4>=.7  \  C<6)*.6  \  GOSUB  2290  \  RETURN 

2435  T(4)*100  \  C'.4>=,8  \  C(6)=.6  •.  GOSUB  2290  \  REM  F  \  REM  Cl  TANKER 
2440  T(5)=90  \  C!8)=.065  \  C(9)=.Q77  >.  REM  M 

2445  RETURN 

2455  ON  F  GO  TO  2460.2465. 2470.2475. 2477  \  REM  COE?  AS  FUNCTION  0?  SS  lOC 
2460  T(4)=100  \  T ( 5)=?5  \  Ct3>=.025  •  C'.9)=.ll  \  RETURN  \  REM  STERN 

2465  T(4)=100  \  T(5»=80  \  C(8;=.029  \  €• 9)=.0?t  \  RETURN  \  REM  SS  AFT 

2470  T ( 4 ) =90  \  T(57=90  \  C(8)=.04  «  C<?>=.07  \  RETURN  \  REM  ON  MS 

2475  T<4)*80  \  T(5)*105  \  C<8>=.054  \  C(9)=.034  \  RETURN  V  REM  FORWARD  MS 
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2477  T(4)=75  \  T(5hll0  \  C(8)*.12  \  C(9)*.0i8  \  RETURN  \  REN  BON  SS 

2480  REN  TANKER  SUBROUTINE  TO  SET  ASSUMPTION  GROUP  44464464644464444466464 

2485  IF  MIDa'Y*  THEN  K=6 

2490  IF  M(2)=,N*  THEN  K=1 

2495  RETURN 

2500  REN  CALCULATION  INITILIZATIQN 

2510  REN  SECTION  2  PRELIMINARY  CALC'S  PRIOR  TO  LOOPS 

2520  REN  CURRENT  41464646464114466 

2530  rem  fx  mmmtntmmm 

2540  C(0»M  \  y  1=1 . 40000E-05  \  R0=1 . 9876 
2550  A- .54 ( 1 • 68354V ( 0 ) ) '  2SR0  \  REN  1/2  RO  V'2 
2560  F(1)=A6A(1)6C!0)  \  REN  PROP 
2570  Sl=1.74Tl»Llt354Dl/Tl 
2580  F(3)*A6S1  \  REM  FRIC 

2590  R1  *1 . 68354*7 < 0 ) 4L 1  /VI  \  REM  REYNOLDS  U/0  F  (THETA) 

2600  F(5)=A6T16B1  \  REN  FORM 
2610  REM  FY  4444444444444144 
2620  C(2)=1+1/((H1/T1)'2-1! 

2630  F(9)*A4L14T14Cl2) 

2640  REN  N  44444444444444 
2650  H<1)=F(9)4L1  \  RESTORE 
2660  FOR  D=0  TO  1 

2670  FOR  C=0  TO  4  \  READ  B(C»D)  \  NEXT  C  \  NEXT  D 

2680  REN  HAVE  DRIFT  44444444444(44444444444444444444444 

2690  REN  READ  IN  TRANSFER  FUNCTIONS 

2700  FOR  D=0  TO  2  \  N‘5  \  P.EM  0= FX  1=FY  2*N 

2710  FOR  C=2  TO  N 

2720  READ  «(C,D)  \  ICXT  C 

2730  FOR  C=1  TO  N 

2740  READ  P(C»D)  \  NEXT  C 

2750  y(l,D)=1.00000€-03 

2760  REN  EXTEND  END  POINT 

2770  K(N+lfD)«12 

2780  P(N+l»D)*PiN-D) 

2790  NEXT  D 

2800  REN  CONVERT  LANHA/L  VALUES  TO  U 
2810  FOR  D-0  TO  2  STEP  2 
2820  FOR  C=2  TO  N 

2830  H(C»D)=SQR(32.2424fI/U.«^;»D)? 

2840  NEXT  C  \  NEXT  D 

2850  REN  CONVERT  M  VALUES  TO  V 

2860  FOR  C=2  TO  N 

2970  U(C»l)sSQR(32.24UtCf:>/Tl’ 

2880  NEXT  C 

2890  REN  NON  ALL  TRANS  FUNCTION  ARE  IN  THE  U  DOMAIN 
2900  REN  CALCULATE  DRIFT  FORCES  M/0  F ^ DIRECT ION ) «**»**« 

2910  C*32.24R04L1 

2920  REN  FX  DRIFT444444444444444444* 

2930  N*0  \  GOSUB  650  \  KM  SPEC  X  TRANS 
2940  FU0)=C4B 

2950  REN  FY  DRIFT  444(4*44444444444444444 
2960  H=1  \  60SUB  650  \  REN  SPEC  X  TRANS 
2970  F(12)=C4B 

2980  REN  N  DRIFT  44444444461444(44144144)4444 
2990  N=2  \  GOSUB  650 
3000  N(2)SC6B6L1 

3010  REN  CALCULATE  WIND  FORCES  U/C  F C DIRECTION ) ************* 

3020  V<4)»(33/H(1»"U/7>«V(1> 

3030  R2*2.40000£-03 
3040  A».54R24(V(4)41.6835'',2 
3050  REN  FX  UIND  4  66  4 
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3040  F(H)*A*A<5> 

3070  REN  FT  WIND  «  t<  t 
3080  V(2)*7/8t<H(2)/33r(l/7) 

3090  V{3)s7/8t(H(3)''(8/7)-H(2}'(8/7))/(33''<l/7))/<H<3)-H<2)) 

3100  C(5)=.92i<V<2)^IA<2)+V(3r2*A(3>)/Ut(2)+R<3)) 

3110  F(16)=A*<A<2)+A<3))*C<5> 

3120  REN  N  MIND  ttWMttmmmt 
3130  N<4)=A4(A(2)+A(3))IL1 
3140  RETURN 

5180  ERASE-TEXT ('TEXT')  \  NO VE. CURSOR* 1.1)  \  RETURN 

7230  PRINT  \  PRINT  \  PRINT  \  PRINT  ..‘FREES  RETURN  TO  CONTINUES  \  INPUT  04(7)  \  RETURN 
32000  END 
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MOOT,  tabular  output  program. 

10  REH  PROGRAM  'STATH00R*  STATIC  MOORING  ANALYSIS 
IS  REH  BY  JAMES  V.  COX  CODE  L44  CEL 

20  REH  26  AUGUST  1981 

21  ren  output  progran  mmmramummtimiitmmmit 

22  REN 

23  REN 

24  REM 

25  COMMON  F(19),T(6)»C(11),H(5)»U(3),D<3)»A(5),V(4),H(3) 

30  COMMON  N(6,2),P<6,2),B(4,l).Dt(7)»Bt(D 

35  COMMON  Rl,Pl,Tl,DlrLl,Sl,Nl,Bl,H0,E0,£l,E2,M0,Hl,T0 
40  GO  TO  8200 

5180  £RAS£-T£XT('TEXT')  \  MOVE-CURSOR! 1.1)  \  RETURN 

7230  PRINT  \  PRINT  V  PRINT  .‘PRESS  RETURN  TO  CONTINUE.*!  \  INPUT  Bt<7)  \  RETURN 

7250  ism  force  output  table  tttmtttmttttmmttntMtmtMmMtt 

7252  Ats*'LLLLLLL  'LLLLLLLL  tttt.ttt.  tttt.ttt.  tttt.ttt, ttt.  tttt.tt* 

7255  PRINT  tB, 'ENVIRONMENTAL  FORCES  RELATIVE  TO  VESSEL  AT  EQUILIBRIUM’ 

7257  PRINT  #B»‘ -  LOCAL* 

7260  PRINT  #B> *  FX  LBS  FY  LBS  M  FT-LBS  ANGLE  (DEGREES)* 

7270  PRINT  #8. 

7280  PRINT  tB.USING  At, ’CURRENT*. ’FRICTION* »F<4) 

7290  PRINT  #0. 

7300  PRINT  *B. USING  At,*  *.*F0RM*.F(6)»F(8).M(0) 

7310  PRINT  *0. 

7320  PRINT  tB.USING  At.’  ■.*PR0P*.F(2) 

7330  PRINT  tO. 

7333  PRINT  tB.USING  At.*  ■.*T0TAL*.F(7).F(8).M(0).T(1) 

7335  PRINT  to. 

7340  PRINT  tB.USING  At.*WAVES*.*<N£AN)*.F(ll>.F(13>»H(3>»T(2) 

7350  PRINT  tO. 

7360  PRINT  tB.USING  At, 'HIND*.*  *,FU5),FU7),N<5>,T(3) 

7370  PRINT  tB,  \  RETURN 

7380  rem  coef  output  table  mtnmmnntmttmtnmnmmmmmm 
7385  At=* 'LLLLLLL  'LLLLLLLL  t.ttt  t.ttt  ttt.  t.ttt'L  t.ttt'L  t.ttt  ##*.* 

7390  PRINT  tB, ’COEFFICIENT  HAXIMUMS’ 

7395  PRINT  tB,’ - * 

7400  PRINT  tB.*  CXB  CXS  THETA2  CY  CMB  CMS  THETAM* 

7410  PRINT  tO, 

7420  PRINT  tB.USING  At, *UIN0V  ’,C(4),C(6),T(4),C(5),’  *»C(8),*  ',Cl9).T(5> 

7430  PRINT  tO, 

7440  PRINT  tB.USING  At, ‘CURRENT*, ’FRICTION* ,C( 1 > »C< 1 ) 

7450  PRINT  tO, 

7460  PRINT  tB.USING  At.’  *, ’FORM*, .1,.1,0,C(3)»**’, CU1), *t’ 

7470  PRINT  tO, 

7480  PRINT  tB.USING  At,’  ’ » ’PROP’ »C(0) ,C<0) 

7490  PRINT  tB, 

7500  At=*t  -  INCLUDES  ANGLE  FUNCTION.  DEPTR/DRAFT  FACTOR  =  t.ttt’ 

7505  PRINT  tB.USING  At,C(2) 

7510  PRINT  tO,  \  RETURN 

7520  ren  results  output  tmummtmmmmmmtmmmmmmmtt 

7530  PRINT  tB, 'STATIC  DESIGN  RESULTS* 

7535  PRINT  tB,* - * 

7540  PRINT  tO, 

7550  PRINT  tB.’HANSER  LOAD  HAUSER  ANGLE  SHIP  ANGLE' 

7560  At=*tttt,ttt.  LBS  tttt.t  DEGREES  tttt.t  DEGREES' 

7565  PRINT  tB,USIN6  At,F(0),T(6),E2 
7570  RETURN 

8000  rem  vidio  output  cgntrol  routine  ttmmtmtmtmutmttmmmttm 

8010  GOSUB  5180  \  REM  CLEARS  SCREEN 
8020  PRINT  „ 'VIDEO  OUTPUT* 

8025  B«0  \  REM  TERMINAL  DEVICE  NUMBER 

8030  PRINT  LI  \  PRINT  8t<0>*  *Bt<D  \  PRINT  Lt  V  PRINT 
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8044  60SUB  7380  \  REM  CQEF 

8054  GOStlB  7230 

8055  SOSUB  5180 

9060  SOSUB  7250  \  REN  FORCE 

8070  SOSUB  7230 

8075  SOSUB  5180 

8080  SOSUB  7520  \  REN  RESULTS 

8090  SOSUB  7230 

8095  RETURN 

9200  ren  output  from  master  MENU  tmumimmtmtmmimttmmmim 
8210  SOSUB  5180 

8220  PRINT  'VIDEO  OR  PRINTER  OUTPUT  (V/P)'i  \  INPUT  D*(7) 

8223  LB=* _ ' 

8230  IF  DS(7)='V'  THEN  SOSUB  8000 

8240  IF  DB(7)=*P'  THEN  SOSUB  9000 

8250  CHAIN  'HOGfi.BAC'  LINE  9000 

8260  REN  VESSEL  OUTPUT  mMMMMMtmmn 

8270  PRINT  #Bf* VESSEL  INPUT1 

8280  PRINT  #B>  * - 1 

8290  PRINT  #B » *  VESSEL  GROUP!  *D«(0> 

8300  IF  DU)=4  THEN  SOSUB  8400  \  REH  CARGO  SUBJECTIVE  INPUT 
8310  IF  D(l)s5  THEN  SOSUB  8450  \  REN  TANKER  SUBJECTIVE  INPUT 
8315  PRINT  BB. 

8320  PRINT  BB. 'DISPLACEMENT  HATERLINE  LENGTH  BEAN  DRAFT' 

8325  AB^'BBB.BBB  LONS  TONS  BBBB.BB  FT  ***.4*  FT  BB.BB  FT' 

8330  PRINT  BB.USING  AB.D1.L1.B1.T1 
6340  PRINT  BB, 

8350  PRINT  BB.'HIND  AREA  INPUT  (AREAS  IN  SQUARE  FEET)' 

8355  PRINT  BB.'END  HULL  SIDE  SUPERSTRUCTURE  SIDE' 

8360  A«='BB,BtB  «*.***  #*»##♦• 

8365  PRINT  BB.USIN6  A*,A!5);A(2),A(3) 

8370  PRINT  BB.'NAX  HEIGHT  ABOVE  HULL  SUPERSTRUCTURE' 

8380  AB='HATER  LINE  BBB.BB  FT  BBB.BB  FT* 

8385  PRINT  BB.USING  A».H(2),H<3) 

8390  RETURN 
8400  REH  CARGO 

8405  PRINT  IB. 'SUPERSTRUCTURE  TYPE  LOCATION  CLUTTERED  BECKS’ 

8410  AB-*  'LL  ' LLLCLLLL LLLLLLLLLLLLLLLLt LL L LLLLLL  V 

8420  PRINT  BB.USING  A<,M(3).M(4).DB(5) 

8430  RETURN 
8450  REN  TANKER 

8460  PRINT  BB.’CI  TANKER  FULLY  BALLASTED  SUPERSTRUCTURE  LOCATION  CLUTTERED  DEI 

9470  AB*‘ 'L  'L  'LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL  'L' 

8480  PRINT  BB.USING  A*,DB(l),D»‘2),DBl4),DB!5) 

8490  RETURN 

8600  REN  ENVIRONMENTAL  OUTPUT 
8610  PRINT  BB. 'ENVIRONMENTAL  INPUT’ 

8620  PRINT  BB.' . 

8630  PRINT  BB. 'CURRENT  SPEED  INCIDENT  ANGLE  HATER  DEPTH1 

8640  AB-'  BB.BB  KNOTS  BBB.BB  DEGREES  BBBB.  FT' 

8650  PRINT  BB.USING  A«.V(0).U<t!.Hl 
8660  PRINT  BB. 

8670  PRINT  BB.'HAVE  SIGNIFICANT  HEIGHT  PERIOD  KNOVN  SIGNIFICANT  PERIOD  INCIDENT  ANGLE* 
8680  Al=‘  BB.BB  FT  'LLL  BBB.BB  SEC.  BBB.BB  DEGREES' 

8685  PRINT  BB.USING  A».H0.DB(6).T0.U<2) 

8690  PRINT  BB. 

8700  PRINT  BB.'HIND  SPEED  INCIDENT  ANGLE' 

8710  AB*'  BBB.BB  KNOTS  AT  BB.B  FT  ABOVE  THE  HATER  BBB.BB  DEGREES’ 

8720  PRINT  BB.USING  AB,V(1).H(1).U(3) 

8730  RETURN 

9000  REN  PRINTER  OUTPUT  CONTROL  ROUTINE  tmtmtMMMMWMtWmm 


SUPERSTRUCTURE' 
BBB.BB  FT' 


CLUTTERED  DECKS’ 


CLUTTERED  DECK* 


I* 


9005  OP®  'LP!‘  fOR  OUTPUT  AS  FILE  t 

9007  8=1  \  REA  PRINTER  DEVICE  NUMBER 

9010  PRINT  #1 t 'PROGRAM  STATNOOR  -  STATIC  MOORING  ANALYSIS' 

9030  PRINT  fl*DAT$'  *CLKi 

9040  PRINT  *1. 

9053  PRINT  U.BilO)'  'B»(ll 

9054  PRINT  Ol.L* 

9040  PRINT  #liTAB(35)» 'USER  INPUT' 

9070  PRINT  IlfLt 

9060  GOSUB  8240  \  REN  VESSEL  INPUT 

9085  PRINT  »1, 

9090  GOSUB  8400  \  REN  ENVIRONMENTAL  INPUT 

9095  PRINT  OlfLi 

9100  PRINT  »lfTAB(35)!'0UTPUT' 

9110  PRINT  *1.L* 

9120  GOSUB  7380  \  REN  COEF 

9130  PRINT  tlf 

9140  GOSUB  7250  \  REN  FORCE 

9150  PRINT  *1. 

9140  GOSUB  7520  \  REN  FORCE 

9200  RETURN 
32000  END 
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MGRAPH,  graphic  display  program 


10  REM  PROGRAM  'ST AT MOOR’  STATIC  MOORING  ANALYSIS 
15  REM  BY  JAMES  V.  COX  CODE  L44  CEL 

20  REM  26  AUGUST  1981 

21  rem  graphics  output  program  mtmmmmmmttmmm 

22  REM 

23  REM 

24  REN 

25  COMMON  F(19)iT(6)»C(11)iM(5)»U(3)»D(3)»A(5)iV(4)»H(3) 

30  COMMON  yi6.2),P(6,2>,B!4,l>.B4<7),M(t) 

35  COMMON  Rl,PlrTl,Dl,LlfSl,Hl.Bl»H0,E0,El,E2»H0,HM0 
10000  DIM  X(19bY(19),X2(14),Y2<14) 

10010  DEF  FNR(X)*XtPI/180 
10020  REM  DEFINE  SHAPES  ItMtttttt 
10030  REM  SHIP-HANSER 

10040  FOR  N=1  TO  6  \  X(N)=2.5*(N-1>  \  X(15-N)=-X(H) 

10050  Y(N)=-X!N)"2/6  \  Y(15-N)^Y(N)  \  NEXT  N 
10060  X<7)=X(6>  \  X<8)=X(9 )  \  Y(7)--100  \  Y(8)=Y<7) 

10070  X(0)-0  \  Y (0) =60 
10073  X( 18)=0  \  Y ( 18) =20 
10076  X( 19 )=0  \  Y(19)=40 
10080  REM  ARROW 
10170  REM  LABLE  COORDINATES 
10180  FOR  N=15  TO  17  \  X(N)=0  \  NEXT  N 
10190  Y<  15)=110  \  Y(  167=123  \  Y<17>--136 
10220  REM  ROTATE  SHIP 
10230  T=FNR(E2) 

10240  FOR  N=1  TO  14  \  SOSUB  10600  \  NEXT  N 
10250  REM  GRAPH  SHIP  AND  HAWSER 
10260  WINDOW!  ,-250, *160. 250, 160,0) 

10265  REGION! 'FULL', 2) 

10270  GRAPH<,-GRIDfLINES,r7»X2U>rY2(lM) 

10275  GWWH,-GRHM.INES,,7,X2<8>,Y2($>,,(2> 

10277  T«FNR!T<6>) 

10283  N*18  \  SOSUB  10650  \  POINT!, X(Y) 

10287  N*19  \  SOSUB  10650  \  POINT!. X,Y) 

10288  N*0  \  GOS4J8  10600  \  POINT! ,X2(0),Y2!0» 

10290  PWBtO)*’  *698(1) 

10300  LMEL(«P1) 

10330  REM  AID  LAHES  TO  ARROWS 
10335  IF  V(0>«0  THEN  80  TO  10365 
10340  N*15  \  T«FMt(U(l))  \  80SUI  10650 
10350  NAP.T!UEXT<X,Y»X,Y> 

10360  MTEXT(,X»Y,*C') 

10365  IF  HO-O  THEN  90  TO  10385 
10370  N»16  \  T«Flt(0(2))  \  OOSUt  10650 
10375  NAP.TO.TEXT (X,Y»X» Y) 

10300  HTEXTlfXfYf'O') 

10305  IF  V(1>«0  THEM  60  TO  10420 
10390  Nil?  \  T4HKB(3>)  \  60SU8  10650 
10400  MAP-TO-TEXT(X,Y?X,Y) 

10410  HTEXTt»X»Y,*«*> 

10420  REM  ANGLE  LMLES 
10430  NAP.T0.YEXT(0,155tX,Y> 

10440  HTEXT(»X,Y,'0*) 

10450  MAP-T0.TEXT(-160,0.X,Y) 

10460  KTEXT!»X»r»*90') 

10470  HAP.T0.TEXT!0,-l55.X,Yi 
10475  HTEXTWX.Yf'lSO') 

10480  MAP.T0.TEXT!155,0,X.Y) 

10490  HTEXT <  ? XtY ?  *270* ) 

10500  FOR  N*1  TO  5  \  PRINT  \  NEXT  N 
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10505  IF  M(7)='DEH0'  THEM  60  TO  10675 

10510  PRINT  'PRESS  RETURN  TO  CONTINUE. *i  \  INPUT  M(7> 

10515  ERASE_6RAPH(»»»0) 

10520  CHAIN  'HOOR.BAC'  LINE  9000 
10600  REM  ROTATION  ROUTINE 
10610  X2(N)=X<N>*C0S(TW(N)I-SIN(T) 

10620  Y2(N)*X(N)*SIN<TW(N)*C0Sm 
10630  RETURN 

10650  XsX(N)*COS(T)+Y(N)t-SIN(T) 

10660  Y*X(N)IS1N(T)TY(N)IC0S(T) 

10670  RETURN 
10675  REN  DEHO  OPTION 

10680  PRINT  ‘NOULB  YOU  LIKE  TO  RUN  ANOTHER  EXAMPLE  (Y/N1S  \  INPUT  D*<7) 

10690  IF  D»(7)=’N'  THEN  GO  TO  10515 

10700  W(7)='DEM0'  \  GOSUB  10800  \  GO  TO  10220 

10800  REH  INPUT 

10805  PRINT  'INCIDENT  CURRENT  ANGLE  IN  DEGREES  *  *5  \  INPUT  UU) 

10810  PRINT  'INCIDENT  NAVE  ANGLE  IN  DEGREES  =  *5  \  INPUT  U<2> 

10820  PRINT  'INCIDENT  HIND  ANGLE  IN  DEGREES  =  S  \  INPUT  U<3) 

10830  PRINT  'VESSEL  ANGLE  IN  DEGREES  =  *f  \  INPUT  E2 
10840  PRINT  'HAUSER  ANGLE  IN  DEGREES  =  *1  \  INPUT  T(6> 

10845  ERASE.TEXT 
10850  RETURN 
32000  END 
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Great  Lakes.  IL.  (  ode  120,  Oakland  CA;  Code  120C,  (Library)  San  Diego.  CA;  Code  128.  Guam;  Code 
154.  Great  Lakes,  IL;  Code  200,  Great  Lakes  IL;  CO  (Code  613).  San  Diego.  CA;  Code  400.  Great  Lakes, 
IL;  Code  400,  Pearl  Harbor.  HI;  Code  400,  San  Diego,  CA;  Code  420,  Great  Lakes,  IL;  Code  420. 

Oakland.  CA;  Code  424,  Norfolk,  VA;  Code  500  Norfolk,  VA;  Code  700,  San  Diego,  CA 
UCT  ONE  OIC,  Norfolk,  VA 

US  DEPT  OF  INTERIOR  Bur  of  Land  Mgmnt  Code  583,  Washington  DC 
US  GEOLOGICAL  SURVEY  Off.  Marine  Geology,  Piteleki.  Reston  VA 
US  NAVAL  FORCES  Korea  (ENJ-P&O) 

USCG  (Smith),  Washington.  DC;  G-EOE-4  (T  Dowd),  Washington,  DC 
USCG  R&D  CENTER  CO  Groton,  CT;  D.  Motherway,  Groton  CT 
USDA  Forest  Service.  San  Dimas,  CA 

USNA  Civil  Engr  Dept  (R.  Erchyl)  Annapolis  MD;  Ocean  Sys.  Eng  Dept  (Dr.  Monney)  Annapolis,  MD; 
ENGRNG  Div,  PWD,  Annapolis  MD 
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